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SYNOPSIS 


In the present work the elastic properties of single 
crystals of Potassium Niohate (crbO^) have been studied hy 
X-ray (thermal) diffuse scattering at room temperature,' In 
this phase, KNbO^ is an orthorhombic perovskite type ferro- 
electric crystal v/ith space group ]Bmm2. The work is pre- 
sented in five sections and four appendices. 

Section one gives a brief account of the theories of 
elasticity, namely the old (Voigt’s) theory based on the 
central force assumptions and the new (Laval’s) theory which 
takes into account the non-central forces in the crystalline 
medium besides the usual central force. The various expe- 
rimental methods available for the determination of the 
elastic constants are discussed, with emphasis on the X-ray 
diffuse scattering method. Then, a summary of the properties 
of KNbO^ is presented, follov.fed by the lattice dynamical 
(Cochran's) theory of ferro electricity of the perovskite 
type crystals. 

Section twn contains a statement of the problem tackled, 
namely, the determination of the elastic properties of KNbO^ 
at room temperature, i.e, in the orthorhombic phase, by. X-ray 
diffuse scattering method and to relate the same to other 
propeirties of KNbO^, 



XV 


Section three gives a simplified formula to obtain 
clastic constants from the diffuse intensity’- obtained by the 
diffractometer and photographic methods. Iho elastic matrices 
for the old theory, having 9 elastic constants, and for the 
now theory, having 15 elastic constants are also given in 
this section. The details of the values of wave normals for 
vaxi^ous directions from a given reciprocal lattice point (relp) 

s 

in the orthorhombic system and their relation to the corres- 
ponding elastic constants are tabulated. 

Section four describes the experimental procedures. 
Specimen preparation is described, followed by a comparison 
of the photographic and diffractometer methods for the dif- 
fuse scattering studies. A simple rotation, oscillation and 
laue camera (-Radius = 5.70 cm), designed and fabrica.ted as 
a part of this work, has been used for tho photographic mea- 
surements. G-. E. XRI)-6 diffractometer with single crystal 
orienter was used for the counter method. Monochromatic 
radiation from Mo target was employed. The various ‘bbrieic- 
tions involved in the methods (photographic and diffracto- 
meter), viz,, general scattering, divergence, skew, polari- 
zation, absorption, second order, mosaic and ?/hite radiation 
and the application of the same to the observed diffuse in- 
tensity have been described. for the estimation of absorp- 
tion correction, the linear absorption coefficient, A » of 
the crystal was determined by the diffractometer method 
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using a scintillation counter to be 55. 9- cm" (Appendix III), 
in excellent agreement v/ith tbe calculated value ( 58.2 
The value of the direct beam intensity was also obtained 
using the Compton scattering from diamond as intermediate 
standard (Appendix II), in order to determine the absolute 
values of the elastic constants. The trigonometrical evsu- 
luation of the thermal wave vectors, the procedure for study- 
ing some properties like equi scattering contours and K-surfaces 
around a given relp, and the choice of the reflecting planes 
for the determination of elastic constants, are also discussed 
in this section. 

The last part i.e. section five describes results obtained 
and a discussion of the same. The values of the (nine) elasi:ic 
constants of FlPbO^, based on the old (Voigt’s) theory are 

= 0.3 (54), 0^2 = 0.5 ( 16 ), C 33 = 0.6 (32), = 3.5 (58), 

C^^ = 3.6 (51), = 2.0 (69), = - 1 * 0 ( 58 ), C^^ = -2.0(27), 

^23 ~ -1.9(14) in units of 10 d 3 mes/cm‘^. Details of these 
calculations are described' in Appendix I. An attempt to obtain 
elastic constants ’■ on the basis of the new theory shows that 
the differences in the elastic constants along certain directions 
obtained from both the theories are not appreciably larger than 
the errors involved in the method. 

The O 33 values obtained by the photographic and dif- 
fractometer techniques employing a very small crystal (volume ru 
1-2 X 10~^ cm^) were found to be nearly identical. Hence 
the time-consuming photographic method has been abandoned for 
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the major portion of the work in favour of the convenient, 
more accurate diffractometer method employing: a proportional 
Counter. Equi scattering contours around the (4,0,0) and (6,0,0) 
nodes are drawn, as also K-sutfaces fot (6,0,0) and (0,0,6) 
relps projected on [OOi] and [OIO] respectively, 0?he relative 
values of the elastic constants obtained from the eq.ui scattering 
contours and IC-surfaces are consistent with the experimentally 
determined elastic constants by X-ray diffuse scattering. A 
discussion of the accuracy of results due to the various sources 
of error is then presented. Ihe accuracy of 0^^, ^225 6^3 is + 10 /q, 

^44» %5’ *^66 ^12’ ^13^ ^23 ~ '^2^, ]?rom 

the values of the elastic constants the Debye temperature of 

KFbO- is estimated to be 302^*K (Appendix IV). A qualitative 
3 

correlation is indicated between thermal expansion and the ther- 
mal diffuse scattering of Potassium Niobate, It shows that the 
directions for which a large spread of diffuse intensity Is 
observed (i.e. for which large amplitude of vibrations of atoms 
exists), are the same as the ones with a large theimal expansion. 



1. JI3TR0DUCTI01T 


The elasticity theories, both the old (Voigt’s) and the 
new (Laval's), are briefly reviewed, followed by a discussion 
of the experimental methods available for the determination of 
the elastic constants. Of these methods. X-ray diffuse scatter- 
ing method is examined in detail. The major characteristics of 
ferroelectric Potassium niobate are summarized. A discussion . 
of the dynamical theory of ferroelectricity is also included. 

1.1 Theory of Elasticity 
1,1.1 Old (Voigt's) theory? 

The usual starting point for elasticity theory is the. 
postulation of Hooke's law, which states that stress T^ is pro- 
portional to strain Ug for sufficiently small strains* ~its gene- 
ralized statement for an anisotropic medium may be taken as 


T. 


6 

E 

S=1 


G H 
rs s 


( 1 . 1 ) 


The constants of proportionality (G^) are called the elastic 

ITS 

constants, stiffness constants or moduli of elasticity. The set 
of linear homogeneous, independent equations represented by 
Bq. (1*1) can be solved for TJ's in terms of the T’s, given by • 
the relations 




6 

s: £ 

s-1 




•. ( 1 . 2 ) 
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The components of inverse matrix S„ are called the moduli of 
compliance. The have dimensions of force per imit area, or 
energy per unit volume and e3q)ressed in units, of dynes/cm^, 

p 

and have units of cm /dsmes. 

X; o 

Pirst Cauchy and later Bom (1915) have developed elasti- 
city theory on the basis of an atomistic approach and treating 
interatomic forces as central forces. On the other hand, 

Voigt ( 1910 ) treated the crystalline medium as homogenoous to 
all scales. Both approaches led to the same conclusion, viz, 
the stress and strain tensors are symmetrical. As a conse- 
quence the nine stress and nine strain components would reduce 
to six each. The six by six array of constants would contain 
56 independent quantities in the most general case. This num- 
ber is, however, reduced to 21 by the requirement that the 
matrices be symmetric on interchange of the double indices 

~ condition follows from the existence of the 

strain energy density W. The symmetry of the C with respect 
to interchange of the subscripts is proved by applying the con- 
ditions of compatibility 


6W = 2 T 6U 

r^s 

since 6W is a perfect differential we have 


( 1 . 3 ) 


T, 


rs 


0W 

W 


rs 


( 1 . 4 ) 
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from Eq.s. (1.4) and (1.1) 

dT^ o2,„ 5T 

fi _ r_ oW _ a w _ s _ n 

rs au ” Su tu FOu ~ au “ sr 
s s r r s r 

Assumption of linearity "between stress and strain allows Eq.(l.3) 
to "be integrated directly giving 


W = 



2 

s > 1 




(1.5) 


The num'ber of independent elastic constants will be fur- 
ther reduced by the symmetry operations of the respective crys- 
tal class (e.g. there are only 9 independent constants for the 
orthorhombic classes^ Bhagavantam, 1966). To investigate the 
effect of a particular symmetiy operation on the elastic matrix 
one develops expressions for the strains in a transformed coor- 
dinate system obtained from the original by the symmetry opera- 
tion. The expression for the elastic energy ¥ (Eq. 1.5) in 
terms of the transformed strains is then equated identically 
to the original ¥, The resulting equation between the coeffi- 
cients of corresponding strain products give the relations which 
reduce the nirnber of independent elastic constants. The conven- 
tional formulation of classical elasticity suffers from the 
disadvantage that the strains are not presented in tensor form. 

As a result any transformation of the co-ordinates requires an 
involved treatment. An alternative formulation is due to 
Sokolnikoff (1946) and others. (Wooster, 1938, love, 1944; Nye, 19571 
Krishnan, 1958, Hiitrtington,X'958;' and Bhagavantam, 1966). The-ir 





fomnilation utilizes the ‘ conciseness and economy of the tensor 
representation. 


The expression for strains is 

9u, 


U. 


rs 


X 

2 


( 


dx 


eu 

r ^ s 




) 


( 1 . 6 ) 


To express Hooke’s law in tensor notation it is necessary to 
treat the elastic moduli as the components of tensor of fourth 
order, 


= °rslcL®kl 

The transformation law for such a tensor is in the case of 
simple rotation 


'^rskl ■" *^pqtT ^rp ^sq ^kt ^Iv* 

This expression which represents quite an involved relation 
provides much simpler method for determining the components of 
C’s in new coordinate system. In equations of motion for an 
elastic medium, the forces on an element of volume are given hy 
the divergence of the stress field, 





8T 

OX 


rs _ 5 

7 "^ 



rskl 






( 1 . 8 ) 


where ^ is density of the material. For the particular case of 
an elastic plane wave aone can take Uj^ = exp,i(ci>t - k.x) 
where are an^ilitudes of the vibration component s , Cci) is the 
angular frequency and k is the wave-number vector corresponding 
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to the wavelength, X == 2%/\'k\» The resiilting equations of motion, 
which are called the Christoff el equations, follow 

In this form the equations constitute the basis for ’long 
Wave’ method employed by Born to develop the elastic constants 
from a lattice theory. 

Bbr applications to the actual situations in which the 
elastic constants are determined from plane wave propagation 
(dynamic method), it is usually preferable to transform to a 
co-ordinate system in which the direction of propagation is one 
of the axes, say x^. In this situation all terms in Christoffel 
equation which involve differentiation with respect to coordina- 
tes other than that along the propagation direction drop out. 

The result is 

S v^ 

where are the elastic constants in the transformed co-ordinate 

system and v is the velocity of propagation, Values for v can be 
obtained solving the secular determinant of Bq. (1.9). The elas- 
tic matrix for the old theory is given by Eq, (1,10), 
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^15 

^16 


u2 

■^22 
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^24 
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^26 


V^ 

■^3 
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C35 
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Cvl 
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?(“45^36> 



2vw 

Aji 


°15 

“46 

“35 

|(“45^“36> 

4(0„+055) 
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2wu 
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_h6 

“26 

“45 

■i(h4'‘'“25^ 

■2^“56'''“i4^ 

|’(Ci2+Ggg2_ 

i 

2uv 


( 1.10 


where C are the elastic constants. 

1.1.2 New (Laval's) theory: 

Bom assumes the existence of central forces i.e. the atoms 
inside a crystalline medium are rigid impenetrable spheres. On 
the other hand, La^al (1951 a, b,) considers a crystalline solid 
to consist of positive ions, which are sensibly spherical and 
rigid, distributed in a periodic array in a continuous fluid of 
valence or conduction electrons. Since the positive ions cons- 
titute almost whole of the mass of the crystal, it is the forces 
which act on them that chiefly determine the elastic properties. 
These forces in the first approximation may be considered to be 
constituted of two parts. The first is produced by relative dis- 
placement of the ions and the second is the interaction of the 
ions with the surrounding defomed electronic medium. Only the 
first part may be taken to be sensibly central but not the second 


one 
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Taking these considerations into account laval (1951 a, h) 
developed an atomistic theory of elasticity and concluded that 
the stress and strain tensors will not he symmetrical, i.e. 

T^y ^ Ty^ and U^y f Uy^ (unlike in the old theory Sec. 1.1.1.), 
so in general there will he nine components of each of the above 
stress and strain tensors. Since elastic constants O’ are given 

X o 

hy the relation, Eq., (1.1) 

= 2 Ug where r, s = 1, 2, ... 9 

s 

there will he 81 constants for a most general case, which will he 
reduced to 45 due to the relation c' = G* . Similar results 

xS SX 

were obtained hy Viswanathan (1954) and Raman and Viswanathan 
(1955). It should he mentioned here that the notations used hy 
Raman are different from those of laval (as given hy le Corre, 


1953). 

Both the 

notations 

are 

given 

below 

for 

comparison. 



XX 

YY 

ZZ 

YZ 

ZX 

XY 

ZY 

XZ 

YX 

Laval ’ s 

notation 

1 

2 

3 

4 

5 

6 

7 

8 

9 

Raman’ s 

notation 

1 

2 

3 

4 

6 

8 

5 

7 

9 


We shall here after refer to Laval’s theory as the new theoiy 
of elasticity. 

The matrix of Eq. (1.10) according to old (Voigt’s) 
theory when expressed in terms of the new theory (using Laval’s 
notation) takes the following form. 
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2c;5 


w2 

^23 


°59 

°27 


(C2^+C^7) 

(Ggg+C^g) 


I 

! 

vw 

E 

A31 


°15 

"^76 

°;8 

(C^Q+C^g) 
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! 
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°i 9 

^26 
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^'^ 28 '^'^ 46 ^ 

C ^ 14 +^ 89 ^ 

(O^g+O'g^ 


uv 


( 1.11 


Prom the matrix Eq. (1.11) matrix A can be calculated 

— 1 —1 

and A^^, ^ 22 ^ derived in terms of the new theory of 

elasticity. 


Thus the new theory of elasticity put forward by laval 
(1951-1957) and Viswanathan (1955) based on the principle of 
homogeneous deformation led to 45 elastic constants. Fore re- 
cently there has been intense discussion about the real diffe- 
rence between the old and the new theories, e.g. Laval (1957) 
questioned the soundness of the rotational invariance condition 
which was overlooked by Viswanathan. Le Gorre (1953, 1956, 1958), 
Raman and Viswanathan (1955) and Joel and Y/ooster (1957, 1958) 
have modified the continuum theory in order to reconcile it with 
the lattice theory and suggested the necessity of taking volume 
couples into consideration. The use of unsymmetrical stress and 

strain components leads to the situation that strain energy can 

eioJLCiH^) 

be- altered by rigidly rotating crystal. Krishnan'^ and Rajgopal 
(1961) have pointed out that the presence of volume couples does 
not lead to any departure from classical behaviour. 
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Experimental verification of the differences in elastic 
constants according to the old and new theories of elasticity 
have also run into difficulties. For example, experimental confir- 
mation of deviation from classical theory of elasticity rested 
largely on observation of light diffraction by ultrasonic waves 
in immonium dihydrogen phosphate (AUPja piezoelectric crystal). 
For this crystal Le Corre (1954) claimed a value of 1.36 + 0.11 
for the ratio while Joel and Wooster (i960) derived a 

value of 1.06 + 0.02 (Bergmann's method, 1954) by analysis of an 
elastogram. Jaffe and Smith (I96l) have measured and of 
AEP by the direct piezoelectric resonance and pulse echo methods. 
These reliable measurements indicate that and 0^^ are equal, 
well within the experimental error limits of 0.1%. Further, the 
results reported on ferroelectric Rochelle Salt (Huntington, 1958) 

f*) 

by various workers are ''poor agreement. More recently, Radha 
and Raogopal (1968) observed the 0^^ and of HaBrO^^ crystal 
to be equal within 0.1%, which is less than their experimental 
error. It therefore seems that unequivocal results to verify 
the new theory are not yet available, 

1.2 Experimental Methods 
1.2.1 St at ic me asurement s : 

The first systematic measurements of elastic constants of 
crystals were made by Voigt (l928) using static methods. Most 
of the investigators immediately succeeding him used similar 
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tecliniques. However static metbod is not favoured in recent 
times compared to the dynamic method because of the low accuracy 
in measuring small strains instead of resonant frequencies, 
which can be measured more accurately. Static method givesthe 
isothermal moduli, while the dynamic measurements based on re- 
sonance techniques give the adiabatic moduli. 

1.2.2 Dynamic measurement : 

An excellent summary of the dynamic methods for measure- 
ment of elastic constants is given by Dine (1952). In the sonic 
and low ultrasonic range (below 1 MHz ) they are generally 
characterized by the establishment of a standing wave resonance 
in a system composed, in whole or in part_, of the specimen under 
study. 

High precision can be attained comparatively easily in 
the dynamic methods, since they involve measurements of the re- 
sonance frequency and not the displacements (as in static mea- 
surements). Dor the freely vibrating system the determination 
of the period and the dimensions of the specimen are needed to 
give the acoustic velocity in the direction of propagation, More 
usually, the system is put under forced vibration and the fre- 
quency of maximum response is found. In both kinds of dynamic 
experiment, information on the damping of the mechanical system 
can be obtained almost in the same operation. Dormula for the 
velocity of propagation are less simple and depend on the 
specimen shape. Special correction factors are involved in 



obtaining wave lengths which change in accordance with the mode 
of oscillation and boundary conditions at the ends of the speci- 
men i.e. whether one end is clamped or both are free (Pine, 1952). 
Yibrational studies of plates have also been used to measure 
els.stic parameters, mostly of piezoelectric materials.lhe possibi- 
lities here are more numerous and mathematically more complex. 
Experimentally there is a greater likelihood for mode coupling. 
Por dynamic measurements, as in static methods, big crystals are 
required which are in many cases to be cut in different directions 
with respect to symmetry axes? studies as a function of tempera- 
ture are also very difficult. These are the limitations for its 
successful application. 

1.2,3 Lattice Interaction with Radiation; 

The methods considered so far are limited to the frequency 
range of mechanical vibrations controllable in the laboratory 
_< 1000 Mc/Sec. This is far below the region where dispersion 
begins. Information regarding elastic spectrum can be gained 
from interactions of the specimen with external radiation. Expe- 
riments with X-ray diffraction, neutron diffraction, infrared, 
Raman and Brillcmin scattering are some exac5)les of interaction 
of specimen with radiation. Of these experiments only X-rays 
and -neutrons have been employed to give systematic quantitative 
data on elastic constants or non-dispersive part of the spectrum. 
Neutron diffraction methods also require very large specimen 
(1-2 cm^); this limits the usefulness of this method whenever 
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it is difficult to grow large single crystals or machining them 
for various orientations. 


Coherent diffuse scattering of X-rays, scattered hy inter- 
action with the thermal waves of the lattice, is a very useful 
tool to study lattice dynamics and to give information on the 
atomic force constants in solids. One uniq_ue advantage of this 
method over other methods is that only a single absolute value 
of the elastic constants suffices to estimate all other indepen- 
dent constants on the absolute scale (Prasad and Wooster, 1955 ). 
Por example, in an orthorhombic crystal , there are six independent 
constants according to old theory G 22 ) ^^ 3 ! G 555 *^ 66 ^' 

If we know the absolute value of C^^ by some other method, from 
the fatios (such as ^ 22 '^^ 1 V which are very correct- 

ly known, one 'can find out O 227 etc. as accurately as 
Another advantage of this method is that the sample remains un- 
strained throughout the experiment (Ramachandran and Wooster, 

1951) unlike in the other methods. The theory of the X-ray dif- 
fuse scattering is now described in detail in the following 
section ( 1 , 2 . 4 ). 


1.2.4 X-ray Diffuse Scattering (Theory): 

Excellent summaries of the theory of X-ray diffuse scatter- 
ing are presented by Born (1942-43), James (1950), Slater (1958), 
Wooster (1962) and Cochran (1969) and of the experimental work 
by Laval (1938-1957) >Lonsdale and Smith (1941-42), Sen (1955) 
and Chakraborty (1958). 
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The idea of determination of elastic constants from the 
study of diffuse X-ray reflection originates from the thermal 
theory of Waller (1923, 1925, 1928) and Faxen (1923). According 
to this theory, atoms constituting the crystalline medium vibrate 
continuously at all temperatures, the amplitude of vibration 
increases with an increase in temperature. These vibrations 
resolve into thermal waves which are superimposed on the static 
periodicity of the crystal. Consequently, regular periodicity 
of the crystal is modulated by each thermal wave and gives rise 
to a series of dynamic stratifications corresponding to each set 
of (h k 1) planes. Coherent reflections from these dynamic stra- 
tifications will be possible at angles slightly different from 
the Bragg angle with respect to the static plane. The intensity 
of such a reflection will depend on the amplitude of modulation 
which is given by the amplitude of the corresponding thermal wave. 
The .amplitude of the thermal wave will be smaller if the frequency 
is higher. Thus the intensity of a dynamic reflection correspond- 
ing to a particular thermal wave can be indirectly employed to 
give its frequency whereas the wave length, which is reciprocal 
of the wave vector, depends on the geometry i.e. the angle of 
incidence and reflection (dynamic reflection). The wave length 
and the frequency of a particular thermal wave being known, its 
velocity can be determined. Since the thermal waves correspond- 
ing to small wave vectors can be identified with the acoustic 
waves, their velocity in a particular direction can be related 
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to the elastic constants of the crystal. Now for the above prin- 
ciple to be applied to practical cases for the determination of 
elastic constants, the thermal theory was required in a form 
such that it could explicitly express the intensity of diffuse 
reflections in terms of elastic constants. This was available 
from the elegant mathematical treatment of the problem by laval 
(1943). Essential features of this treatment are described in 
the following paragraphs. 

Thermal vibrations of a crystal consisting of N unit cells 
and g atoms per unit cell can be resolved into 3 Ng independent 
waves. Assuming the thermal waves to be plane, the displacement 
fj. of an atom situated at the extremity of a vector (I + J) can 
be expressed as, 

3Ng 

"d = 2" (V - 5a Cl + 

where , 

I 1^3.^ •••/ 

iT (T^a^ "r ^2^2 «T^a^, (J^, ^3 ^ 

a^, SLg, are the lattice parameters, 
is the frequency of a thermal wave, 
q^ is the wave vector of the thermal wave 
n. is the phase angle at time t flcnd 

A. is the amplitude of vibration. 

3 ^ 

Neglecting the absorption and dispersion and excluding 
the Compton scattering, the total amplitude scattered by a 
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crystal vibrating under the influence of 3Wg tbermel waves is 
given by the expression, 

{ ^ ^ _ 

6 E 2 f . exp. 2u i Ivt + X (l+J) + Z X A. Cos 2 -k 
I D ^ I cc=1 

- 5ci (I+J) +^3011) 

( • -Tl 520 

= C E E f . exp. 2% i i vt + X (I+J) \ 77 exp. 271 i 
13^ (. j a=1 

X Cos 2n S)„t - (I+j) + (1.13) 

where f . is atomic scattering factor, 

^ S - S 

X = ■“ — T ; S and S„ being unit vectors along the directions 

X O ^ 

of reflection and incidence, /\ is the wave length of the incident 
radiationjC is the amplitude of the wave scattered by a free elec- 
tron and is given by Thomson’s formula 

e = (A^ + Cos^ 

mc'^r 

where e and m are respectively the charge and mass of the electron, 
0, diffraction angle, c the velocity of light, r the distance of 
the point of observations from the electron, A and B are the com- 
ponents of incident electric field perpendicular and parallel to 
the plane of diffraction respectively. 

let us put 27 cX. A^ = Z 

and 27t (I+J) + n^ J = 03 
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we then get 

exp. 2n i X Cos 2iz (I + J) + 


Da 

iZ Cos63 


= J^(Z) + iJ.|(Z) (e^"+ ,+ Ag(Z)(e^^^+ ) 

+ + i“ J^(Z)(el”H. 

where, is the Bessel function of the order n. Now, the total 
radiation scattered consists of the following components: 

i) The laue-Bragg diffraction having the same freq_uency as 

that of incident heam and given hy the expression 

g p 2 gi2TtXI 


where , 
and 


B = Z f. H. 

0 2 j 

3Ng 

-"o (Zit X 3- factor which arises 

CX“~ 'I 


due to thermal vibrations and reduces the scatterin,^ factor f. * 

j 

* • \ th 

ii; Laval considered general case of n order radiation due 
to thermal agitation and gave the following expression for the 
total intensity scattered (l.^) in the first order. 





B 


a=1 


where 


and 


J.(2ii:XA. J 
““ D ^ ^ Jq(2ixXI^^) 


-aM' 


gi27i:(MJ + n^.^) 


(1.14) 


B = 2 f . H. . gi2-rt(MJ-. n..„) 

-aL ^ 3 3 J^( 27 tXA.^) 


Da^ 


and ii = (X + q) . This has been represented hy a vector diagram 
( Big. 1.1). 
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Pig. 1.1: Q is tile pole of the diffusion; X reciprocal 
lattice point nearest to it. 


JQ-l other radiations of the first order which do not have 
wave vectors equal to + QX do not produce any diffuse scattering 
in the particular direction under consideration. Por poles very 
near to the corresponding nodes, the fraction of the incident 
intensity scattered in the first order due to wave vector q , is 

Uu 

expressed hy 


(IT = 


e^Njx] 


m 


3Ng 

2 

a=1 


P 


a 




a 


E f . H. 0 . 

V 


Cos IX. A 


da 


(1.15) 


1 1 

where E = — — — — + hy 

hva/kT_ ^ 


a’ 


is the mean energy of 


vibration, m is the mass of the crystalline motif and 0. are 

J ^ 

coefficients proportional to such that, 


A. 


3a 


= i i lg'. = A 

c ! c !* “ 

3a 3a 


A being the mean amplitude, and Em. C. 

^ J J oc 


in 




/ 

/ 
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For polyatomic lattices one wave vector corresponds to 
a numler of thermal, waves having different frequencies. So no 
exact information about the frequencies can be obtained from 
the intensity of diffuse reflections. But for very long waves 
i.e. short wave vectors (q -^0), however only three acoustical 
waves are primarily effective in producing diffuse reflections. 
Since q -» 0 the frequency of the three acoustical waves ( one 
longitudinal and two, transverse) approach zero while that of 
optical branch remains still high. Tbe amplitude scattered by 
optical branch becomes inappreciably lovi/ whereas that due to 
acoustical branch remains high (Fig. 1.2). 



Fig* 1.2: Optical and acoustical branches of the dispersion 
relation for a diatomic linear lattice showing the 
limiting frequencies at q = O' and q = qmax = ti:/2a. 
The lattice constant is 2a and M.|, M 2 -are the 
masses of the two atoms. 

For long waves C. can be approximately taken to be unity, and 

j a 

(n^ - n!’ ) etc. become zero as the amplitude of vibrations 

of different atoms due to wave a can be taken to be same in 
this case. Then the expression (1.15) reduces to, 
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N. . kT 3 




m 


S 

a=l 


Cos^ lX.I 


10. 


or 


Since, E 


a 


kr as la'i) ^ « kT 


t57 


!!izli! IL I is.Aa 


V 


[2 a =1 


V 


( 1 . 16 ) 


a 


where E^, is structure factor of the plane (h k 1) under consi- 


deration 


is density of the crystal, 

■? is volume of the unit cell, 

V ' is the velocity of thermal wave vector q and 
a is the thermal wave. 


As explained above the contribution of the optical branch 
of thermal vibrations to the diffuse intensity is negligible for 
q 0 as compared to the contribution of the acoustic waves, the 
velocities of which can be expressed in terms of the elastic 
constants of the crystal. 


Since, 




Y 2 
a 

11 

A^1 

P“ 

^X 

4 * 

Ai2 

a 

ry 

4 - 

A ^3 

pCC 

2 

V ^ 

a 

ry - 

Ai2 

p 05 

^X 

+ 

^22 

a 

j- 'Y 

4 * 

A23 

P“ 

■ z 

V ^ 
a 

p® _ 

■^13 

a 

^X 

4 * 

■^23 

pOC 

4 * 

A33 

p“ 

^z 


(1.17) 


where, P^, Py, P^ are the polarization vectors of the wave, i.e. 
the components of the imit displacement vector for the wave a. 
The elements A 22 are given by the equation (1.10). 
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Now, 



+ 2 NLP^ Pj 


(I£P“)^ + (NP“)^ + 2 MNPy P^ 

2 !!''[?“ P^ ] (1.18) 


where 1, M, N are the direction cosines of the vector X and 
u, V, w are the direction cosines of a thermal wave vector 
with respect to the orthogonal system of the reference axes. 
Solving for the expressions, 

(IP“)^ IMP“ p“ 

__A— i etc. from Eqs. (1.17) and 

(1.18) we get, 


3 

E 

1 


0os2lx.l 


‘a 5 






= A“!j + A22 + N^ + 2 


a 


+ 2NLA3^ + 2IM“2 


(1.19) 


= K ["u, V, wl 

— 1 — 1 

where A. ! are the elements of the matrix A which is inverse 

X J 

to the matrix 




All 

CM 

Ai3 


A = 

A^2 

^22 

^3 


1 

_Ai3 

^23 

■^33_ 

have 

from ( 1 . 16 ) and ( 1 , 

19) 


(57 

kl 

X^ 

■ . K r 

n XT 

V 


U. V 1 


Ihkl 


(1.20) 
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A similar expression has also been derived by Begbie 
and Born (1947) for q 0 by a q-uantnm mechanical treatment. 


Y/e have till now considered the contribution of 1st order 
diffuse scattering only, laval has shown that the 2nd order 
diffuse reflection from, a pole of diffusion will be produced 
if the resultant of the two interacting wave vectors and q^ 
is equivalent to q as illustrated in the diagram. (Big. 1.3). 



Big. 1.3: Vector diagram for second order 
diffuse reflection. 

The intensity diffusely scattered in the 2nd order is given by 


fT" - y 7 

^ ^ 2T ^ a6 v^ 

^ a 


Cos2 (X.A^) Oos"^ (X.Ig) 


q 


a 


-2 

'ip 


(1.2l) 


Ramachandran and Wooster (1951) made certain approximations in 
view of the fact that CJ^<< 0)pand deduced the following expres- 
sion from Eq. (1.21), 

.3 o V4 Cos^ 


(Ts. = 


n 


R(kT)‘ 

V 


p2 y 

^T 


V^ 


( 1 . 22 ) 


a 


3 Cos^ (X.l^) 


Z 

1 




can be expressed as a function of 


elastic constants just as Eq. (1.19) and we may denote it by 
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k'Hu, v, 

Then one gets. 


= jfl Eimf. It 

2 2f T q 


k'L'j 


y/ 




hhl 


(1.23) 


This expression can he used (if required) for correcting the 
observed intensity for the second order effect before evaluating 
the elastic constants from Eq. (1.20). This point will be dis- 
cussed further in Section 4.5.6. 


The matrix Eq. (1.10) has been expressed in terms of 
classical theory (Voigt's) of elasticity. Now according to new 
theory of elasticity (Sec. 1.2) the matrix Eq. (1.10) taJtces the 
form as expressed by the Eq. (1.11), Erom this Eq. (1.11), the 
matrix A can be calculated and A)^^ , ? ®tc, can be derived 

in terms of the new theory of elasticity and on substitution in 
Eq. (1.19), it can be combined with Eq. (1.10) to express the 
intensity of diffuse reflection corresponding to "particular di- 
rection of propagation of thermal wave and particular reciprocal 
lattice point. The expression in Eq. (1.19) is generally sim- 
plified (involving very small number of elastic constants, even 
one or two in some cases) for certain nodes and simple directions 
of propagation of thermal wave. Observations on intensity of 
diffuse reflections along such directions would yield elastic 
constants. 
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1.2.5 Ihe X-ray diffuse scattering (Experimental): 

The method was for the first time successfully employed 
hy Ramachandran and Wooster (1'9^) for the determination of 
elastic constants of some cubic crystals. They used a Geiger 
counter spectrometer suitably modified for the measurement of 
feeble X-ray intensities. The method was later on used by 
Wooster and his collaborators (Prince and Wooster 1951 » 1953, 
Prasad and Wooster 19535 1956 etc.) for determining the elastic 
constants of a number of crystals (mostly inorganic cubic and 
more symmetrical tetragonal i.e. 4 2, 4 2m, 4/mm - crystal classes). 
They have also determined the ratios of the elastic constants 
(Hoemi and Wooster 1952, Prasad and Wooster 1956) by photogra- 
phic method of measuring intensities. To derive absolute values 
of elastic constants Wooster and collaborators have used bulk 
modulus data taken from other sources (Bridgman., 1938, 1945) 
or Compton scattering from diamond (Compton and Allison 1935). 

Sen (1955), Chakraborty (1958) and Chakraborty and Sen 
(1958) have developed and standardised a method (of the quanti- 
tative study of absolute measurements of the intensity , of diffuse 
reflections by photographic photometry also cietormine the 
elastic constants of Bensi4 (crystal class B^) from X-lray* "study . 
only. 
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1.3 Potassium. Niobate (KHbO^) 

KKTdO^ is an orthorhombic perovskite type crystal at room 
temperature. It exhibits ferroelectricity (Matthias 1949), a 
property of the crystal having spontaneous polarization which 
can be reversed by means of an electric field, ( Ferroelectricity 
was ably summarized by Jona and Shirane (1962), Patuzzo and 
Merz (1967)). 

1.3.1 Phase transitions in EUbO^ crystal: 

A crystal of Potassium Niobate undergoes the following 
phase transitions (Shirane et al(l954) jSchcnk (196 9) ) in the 
same order as in BaliO^. 

Rohombohedral orthorhombic - Tetrago- 
nal cubic. The Ourie temperature of KNbO^ lies at 

435°C. All these phase transitions are of first order and con- 
nected with detectable temperature hysteresis (Shirane et al 
(1954)). 

1.3.2 Single crystal growth* 

Single crystals of KFbO^ were first produced from binary 
melts using an excess of as a flux (Wood, 1951) Molten 

mixture of 1,2 mole of and 1 mole l!rb 20 ^j gives quite 

large flawless crystals (Miller, 1958), The phase relationships 
in the system were studied by Eeisman and HoltzberB 

(1955, 1956). Ihe melting point of KlbO^ is about 1050^0. 
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1.3.5 Dielectric properties: 

The dielectric constaht of multidomain crystals of KlTbO^ 

as a function of temperature were studied hy Wood ( 1 951)? Timofeeva 

et al (1962) and in more detail hy Shirane et al (1954). The 

dielectric constant vs. temperature curve obtained hy Shirane et 

al? shows a marked resemblance to the higher branch of the curve 

for BaTiO, (Merz, 1949). The temperature dependence of the 

spontaneous polarization was measiu?ed by Triebwasser''( 1956) . The 

—6 2 

value- of P_ just below Curie temperature is 26 x 10” C/cm . 

o 

1.3.4 Thermal Expansion: 

The temperature dependence of the lattice parameters of 
KlTbO^ was studied first by Wood (I95l) and in more detail by 
Shirane et al (1954). The general appearance of the curve is 
again quite similar to the corresponding curve of BaTiO^ (Kay 
and Vousden, 1949)? but the spontaneous strain in KHbO^ is lar- 
ger in all the three phases. This is quite expected since 
the Curie temperature which is generally indicative of the 
strength of interactions? is larger in the case of KEbO^ than 
that of BaTiO^. Accordingly transition energies are considerably 
larger for KEbO^ than for BaTiO^ at each of the three phase 
changes. (Shirane et al 1954? Triebwasser et al 1955). 

1.3.5 Piezoelectric and Elastic properties; 

The dielectric anamolies and phase transitions in single 
crystals of Potassium niobate are discussed in sections 1.3.3 
and 1.3.1. We may proceed a step further and enquire which of 
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"the piezoelectric and elastic coefficients "behave anaiaolously 
as the transition temperature is approached? TJn fortunately, for 
this particular crystal no detailed studies have "been made in 
this direction. Dielectric and piezoelectric properties have 
been reported hy Timofeeva and Popova ( 196 2), They have reported 
only one piezoelectric constant d^^ and its temperature depen- 
dence. ho elastic constant data have been reported to date on 
single crystals of Potassium Niohate (as far as author is aware). 
Egerton and Pillon ( 1959) , Jaeger and Egerton (196 2) have studied 
ceramic samples of haQ ^qKq They have reported Young's 

modulus and two piezoelectric constants d^^ and d^^ respectively. 
Their study would have "been useful for extrapolating constants 
for pure KhbO^. Unfortunately they have not reported the data 
for other compositions in detail, Erom theoretical point of view, 
this system (ha, K) FbO^ was studied in detail by Cross (1956), 
who applied the thermod 3 mamic treatment developed by Devonshire 
(1949) for BaTiO^ to both NahbO^ and the system (ha, K)NbO^. 

His theory explains only dielectric properties and the various 
transition temperatures observed in the system. 

1,3.6 ‘Crystal structure: 

The structure of Potassium hiobate at room temperature 
has been investigated in detail by Katz and Megaw (1967). The 
structure is strictly isomorphous with orthorhombic BaTiO^,with 
the space group Bmm2 and 2 molecules per unit cell, but all 
deviations from perovskite aristots^pe are rather larger in 
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KI'TbO^, A good comparison regarding the structure of KETbO^ and 
BaTiO^ in orthorhombic phase is given by the above authors. Cell 
dimensions given by Katz and Fegaw at room temperature are 
a = 5.697 b = 3.971 % c = 5.721 i in agreement within about 
0.001 S. with those of Vousden (1951)? Shirane, Nevmham and 
Pepinsky (1954). Zr filtered MoKa radiation was used in the 
above investigation. Structure factor data used in the present 
investigation are taken from Katz and Megaw (l967) and Dwivedi 
and Srivastava (1969). 

1.3.7 Comparison of KllbO^ and BaTiO^ 

It is worth mentioning that to date KNbO^ is the only 
perovskite ferroelectric crystal that exhibits the same phase 
symmetries and same sequence of transitions as BafiO^, A compa- 
rison between some of the pertinent data of the two crystal, taken 
from Jona and Shirane (1962), are shown in Table 1,1. 

1.4 Lattice dynamical theory of ferroelectricity; 

A very successful attempt has been made by Anderson (1960) 
and in more detail by Cochran (1960?61, ' 67,69) and Cowley (1963) 
to provide a microscopic theory (based on Lick and Overhouser 
model, 1958) of ferroelectricity for certain crystals (Perovskite 
type, displacive ferroelectrics etc). This is based on the assunp- 
tion that ferroelectric transitions are the result of an insta^ 
bility for a particular normal mode of vibration, and they can 
be treated as a problem in lattice dsmamics. Let ails of this 
theory are presented by Patuzzo and Iferz (1967). . 
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Table 1.1 

Coii 5 )arison between KNbO^ and BaTiO 

ZJIbO, 


Transition temperatures 
°C 

Transition energies 
(Cal/Mole) 

Maximum tetragonal 
distortion c/a 

Spontaneous polarization 
at the Gurie point 
(10“^ C/cm^) 

Curie constant C(°K) as 
defined from 6 = o/T-Tq 
Te - lo 


435 , 225 , -10 


190 , 85 , 32 


1.017 


26 

2.4x10' 

58 


Coefficients ^ 


Sl1^ and\^^^x ^ _ 
of the free energy expansion 
2 . 1 >r ■o-n4 


10x10 

c. g.s. 


,-15. 


1 2 1 
A = Ts X XP + 4 
2 2 






?\iii 


XP 


'^^^^X=54x10~^^ 
c. g.s. 


3 

BaTiO^ 

120 , 5 , -90 

49 , 21 , 11 

1.01 

18 

1 . 7x10^ 

11 

-1'^ 

-20x10 c.g.s. 
25x10“^^ c.g.s. 
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Essentially, the theory is based on the argnment that 
if a crystal is wholly or partly ionic, lattice vibrations are 
accompanied by polarization oscillations of equal frequencies 
which create a local field interacting with ions through long- 
range Coulomb forces. If for one particular mode of vibration, 
these long range Coulomb forces and short range forces balance 
out, then the crystal becomes unstable for this mode. 

The detailed mathematical treatment is not given here. 
Only the final results, as a consequences of the theory, are 
quoted below. 


The dielectric constant C is connected to the frequency 

s 

of a critical mode ty Bl, (1.24) 

(1.24) 

where K is a proportionality constant and is independent of 
temperature. As the temperature approaches the Curie point, 
becomes extremely large and tends to decrease rapidly. The 
occurrence of ferroelectricity and the crystal instability can 
be connected by lyddane, Sachs, Teller (I.S.T.) relation (1941)- 



(1.25) 


€ the static dielectric constant, 
s 

and optical dielectric constant is equal to the 

square of the refractive index. is longitudinal frequency 

andCJgi the transverse optical frequency at q = 0 (wave vector). 
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This L.S.T. relation is written for two atoms per nnit cell. 
Cochran (i960) has extended this treatment for a more general 
case of n atoms in the elementary cell. 


e /€ 


e 



( 60 -j ) gn 


( 1 . 26 ) 


(CO ) and ( tO.)mha-ve the usual meaning as stated above. This 

J JJ J 4_ 

relation is obtained from lattice dynamical equation (Cochran 
i960). The equation (1.26) is derived on maiy ph3i'sicclly 
unrealistic assumptions; one of them is that this relation is 
true only for diagonally cubic crystals i.e. ones in which every 
atom has surroundings of the tetrahedral symmetry. However, no 
ferroelectric crystal having this symmetry has yet been disco- 
vered (Hatuzzo and Merz 1967). Further shell model treati^ent is 
not adequate in the electronic frequency range (though these fte- 
quencievS do not appear in the final results). Further we assume 
that the Eq. (1.26) (although this has not been proved conclusively) 
is valid for every cubic crystal. If one of the transverse opti- 
cal branches vanishes in the neighbourhood of q = 0, at a certain 
temperature, then the static dielectric constant (Bq. 1.26) be- 
comes infinite at that temperature which is precisely what happens 
at the transition from ferroelectric to paraelectric state, i.e. 
frequency 6i)p( 2m of a ferroelectric mode varies with 
(T - Tq) as 

O) j a (I - I„) 


(1.27) 
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To is Curie temperature. This equation is generally consi- 
dered a fundamental consequence of Cochran’ s. theory (I 96 O 56 I) and 
one which provides a good check for its validity. 

Since the 'ferroelectric soft mode’ of ahnormsHy low fre- 
quency lies in the far infrared following experiments have been 
used to search for such a mode. 

i, Far infrared (Reflectivity expt, also) 

ii. Raman spectrum 

iii. Dielectric dispersion 

iv, Neutron diffraction 

There are some experimental difficulties which one comes 
across in verifying the existence of soft ferroelectric mode 
(Cochran’s mode) for example, undoubtedly Raman and Infrared 
experiments could clarify the occurrence of low frequency trans- 
verse optical mode (Gochran'-s mode) but because of the limita^- 
tions of the experiments, it is difficult to find the ferroelec- 
tric mode in perovskites by analysing their far infrared data 
(Barker 1966). Details are described by Dvorak (1967). Now 
low frequency mode in millimeter range is experimentally very 
difficult to handle and if a crystal possesses a high conducti- 
vity (as in the case of ENbO^) near the transition temperature 
the dielectric dispersion study becomes impossible. Neutron 
dispersidn ,may be comparatively easy if special care is taken 
for several impoirtant corrections, further this requires very 
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large single crystals. Experimentally the existence of this soft 
phonon mode 'and its strong temperature dependence have been obser- 
ved in SrTiO^ and KTaO^ (Yamada and Ehirane 1969). These are 
both perovskite tjrpe crystals which show an anomalous temperature- 
dependent dielectric constant, but neither of these materials 
actually undergo a ferroelectric transition. On the other hand, . 
in the case of BaTiO^ which uudergccs a ferroelectric transition 
around 125°C, the results on the existence of a soft phonon mode 
as determined by different experimental techniques are found to be 
in disagreement with each other (e.g. Barker, 1966 and Ballantyne, 
1962). Recently, Harada et al (1967) deduced a soft phonon mode 
energy from X-ray diffuse scattering studies. The significance of 
these X-ray results is the fact that this optic mode energy, lies 
much below the energy of the acoustic mode. This ba.s been confir- 
med by the thermal neutron study of Shirane et al (1969). It is 
interesting to note that Bell and Ruppercht (1963) have observed a. 
minimum in the elastic constant vs. temperature curve of SrTiO^, 
but no corresponding maximum or even change of the slope in the 
dielectric constant vs. temperature curve. If the particular 
short range lattice forces that determine the elastic properties 
are also connected with the frequency of the soft mode, a disconti- 
nuity of dielectric constant would be expected. These results 
how that the short range forces responsible in SrTiO^ for sound 
opagation are not connected with dielectric properties. It is 
surprising to see that a marked change in the elastic 
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constant ys. temperature graph can occur vd.tli little or no obser- 
vable effect on the dielectric constant. Similar observations 
have been reported by Cummins (1970) on Gd 2 (MoO^)j. 

Recent worh due to Shapiro and Cummins (1968) on phase 
transition in Quartz, imoros (1968‘) on lTa]\f02 and triglycine sul- 
phate, by Hon jo and Harada (196 9) on a, and p, Qurrtz, shows that 
optical mode of low freq.uency is not is sponsible for the ferro- 
electric transitions. (Cross et al 1968) (T,G. S. and Quartz are 
the examples cited by Cochran (1960-61) in his theory of displa- 
sive Ion model). Very recently, a disorder theory for BaTiO^ and 
OlbO^ is proposed by Lambert, Comes and G-uinier ( 
in which they claim to account for the ferroelectric transition 
in a better manner than the alternative explanation of soft mode by 
Do chran ( 1 96 0-6 1 ) . 

The theory of ferroelectricity developed by Cochran has 
been briefly summarized above. How we shall discuss how the elas- 
tic constants of KFbO^ determined from the X-ray diffuse scatter- 
ing studies are likely to be a.ffected in the light of Cochran’s 
theory. 

If one recalls the expression for the diffuse i-ntensity, 
it is comprised of two parts, intensity of acoustic vibrations 
and the intensity due to the optical vibrations. 


^diffuse ^ac ^opt 


(1.28) 
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Generally the optical contrihution is negligihly small (owing to 
high frequency) and hence is eliminated as general scattering 
correction. If the elastic constants of this crystal (KETbO^) are 
known hy methods other than X-ray diffuse scattering, one can 
easily find out the contrihution of the low lying optical fre- 
quency and from this the corrected elastic constants. But to date 
no such elastic constants data are available (as far as author is 
aware). Secondly if the infrared or Raman spectrum of this com- 
pound is available then one can try to find out the optical fre- 
quency contribution to the total diff^ise intensity. Recently 
perry et al (1969) have reported the Raman laser spectrum of 
Potassium Niobate, in all the four phases. They also studied 
the temperature variation in the rhombohedral and orthorhombic 
phases. But in their spectrum no such low frequency mode has 
been observed. The temperature dependence of phonon associated 
with different frequencies is hardly 1%. The lowest frequency 
in the orthorhombic phase which they have reported is 298 cm . 
Potassium Mobate has a high Curie temperature, around 708°K. 

Then if we recall the famous Cochran's equation (1,27) we see 
the difference T-Tq is very large at room temperature and 
would be definitely finite and hence the elastic constants will 
not diverge from their normal value. If one makes a very crude 
approximation that KHbO^ is cubic at room tem-perature and finds 
the corrections to the elastic constants as shown by Cowley (1964), 
the correction does not exceed normally 1-3% . Bae account 
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cpjn be talvcn of euch correctioias for lovv frequency opticoJ. 
mode, Thai: is, the values of the elastic constojits of a 
single crystal of Potassium Niobate would be hardly affec- 
ted <at 1*00111 temperature due to this trojis verse phonon mode 
(Cochran's mode) of low frequency. 
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2. STATEl/ENT OP THE PROBLEM 


The behaviour of a ferroelectric crystal is characterized 
by its dielectric, elastic and piezoelectric constants* Such 
complete data are available on very few ferroelectric materials. 
Potassium Niobate with its high Curie temperature (708°K) is a 
potentially useful ferroelectric. The interest in this compound 
is further enhanced by its close similarity to. the much studied 
ferroelectric Barium Titanate (Table 1.1). However as pointed 
out in Sec, 1.3.3, fairly detailed dielectric data are available, 
while Icnowledge of the elastic and piezoelectric constants of 
Potassium Hiobate is lacking. The main purpose of this investi- 
gation therefore is to evaluate the elastic constants of Potas- 
sium Hiobate (EPTbO^) at room temperature i.e. in the orthorhombic 
phase. X-ray diffuse scattering method has been employed for this 
purpose in the present study, both by photographic and diffracto- 
meter technique. 


A crystal of Potansium Hiobate belongs to an orthorhombic 
\im2^ class'^ has 9 elastic constants according to old (Voigt’s) 
theory and 15 elastic constants according to new (Laval’s) theory. 
In order to measure all the nine elastic constants (old theory) 
a crystal under investigation has to be mounted along all the 
three mutually perpendicular axes. The atudy of the following 
Bragg reflections were made on the XRB-6 (G.E,) counter diffracto- 
meter with CA-7 MoKa target, operated at 35 KVP/25 mA, 
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Relp (Reciprocal lattice Point) (6,0,0) witli c-axis verti- 
cal (i.e. perpendicular to the equatorial plane of the diffracto- 
meter.) gives the constants and independently and the 
dependence of constant other constants. 

The relp (0,0,6) with h-axis along the vertical gives the 
constants 0^^ and 0^^ independently and the dependence of 0^^ on 
other constants (Tables 5.1 , 3.2 and 3.3). 

The relp (0,4,0) with a-axis vertical gives the constants 
Cgg and independently and the dependence of on other 
constants. 

The choice of the above planes (relps) for the study is 
explained in Section 4.5.10. 

In this way all the constants (of the old theory) have 
been obtained (Tables 3.1, 3.2, 3.3). Additional constants of 
the new theory also have been evaluated and shown to be similar 
to the ones obtained from the old (Yoigt's) theory. In each relp 
normally three rekhas (any line drawn through the given relp) 
have been studied on either side of the relp. For each rekha 
5 to 7 sets of observation are recorded on the proportional 
(Xe-filled) counter. For each observation, time for 10,000 
counts was recorded. Zr-Y balanced pair of filters have been 
used during the entire diffractometer observations. On the 
other hand photographic data were obtained on rotation, oscilla- 
tion and laue camera with crystal (liF) reflected monochromatic 
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i'-oKa radiation (operated at 35 KYP/25niA). Each exposure for 
diffuse spots, ran for 35 to 40 hours. Type F X~ray films 
were employed. 

In addition to this the relps (4,0,0) and (0,0,4) with 
h-axis vertical have been studied on the XKD-6 diffractometer. 

The crystal used for these relps was very small (Tol, 1.782 x 
10”"^ cc). The res^iLts obtained from these relps are compared 
with those obtained by photographic technique on an equally 
small crystal (Vol. 0.996 x lO”'^ cc) from the relp (4,0,0) with 
b-axis vertical. The constants obtained by these two different 
techniques are in very good agreement, suggesting that one can 
also work on the diffractometer with such a small crystal, as 
the one used in the present investigation. The diffractometer 
technique is obviously more convenient and handy, besides giving 
better accuracy. 

Some 'of the thermodynamic parameters can be obtained from 
the observed elastic constants. One such parameter viz, the 
Debye temperature is obtained. An attempt is made to inter- 
relate the elastic constants and thermal expansion as to correlate 
the same with the crystal structure. 



3. A SIMPLIFIED EORMOLl fOR THE DIFFUSE 

INTENSITY AND THE ELASTIC CONSTANTS 


The general formula for the ratio of the intensity of 
diffuse X-ray reflection I^ (1st order only^ second order correc- 
tion is neglected for the present) to that of the incident heam 
I^ corresponding to thermal wave vector q and the reciprocal 
lattice point (h,k,l) is given hy Eq. (1,20). The corresponding 
formula for the diffractometer and photographic methods are given 
hy Bq. 3.1 and 3.2 respectively. 


3.1 Foimulae 


Diffractometer method ( Ramachandran and Wooster, 1951) 






2 ip ^ C 

T -2 1- 

per unit cell. 


hkl 


(3.1) 


where 6 


„2 


(A^ + Cos^ 0) 

m. c‘^ 


Symbols have been explained in Sec, 1.2.1. 


Photographic method (Srivastava, i960, Joshi 1964? 
Kashyap, 1964) 


rt-"- FT ^2 X^ vrii .. wl 

^1 “ i^ - ;;2 -2 


T 


.2 2 


where e2 = ( — | — ) (a^+B^ Gos^ 0) 

mc'^r 

r is the radius of the camera (r = 5.70 cms.) 


(3.2) 
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Kc is the linear absorption coefficient of the crystal 

6V is the volume of the crystal 

V is the volume of the unit cell 

T is the absolute temperature at which observations 
are made. 

k is the Boltzmann constant 
_n is the solid angle subtended at the crystal by the 
counter collimator slit 
is the direct beam intensity 

is the structure factor of the plane (h,k,l) at 
temperature T. 

1 is the thermal wave vector 

X is the rel vector (a line joining a given point (relp) 
to the origin of the reciprocal lattice) corresponding 
to the given relp (h,k,l) 

and 

K[2u,v,wJ^j^^ = + M^A22 + + 2mA~^ 

+ 2ELA~;[ + 2L¥jql (3.3) ' 


where 1, M, N and u, v, w are the direction cosines of X and 

q respectively corresponding to the orthogonal elastic axes, 
i — i 

A. i are the matrix elements of the matrix A ; 

J- (J 

where, Ia^^ A..^ 


A”" ^ 

det, A 


and det A = A 
A 


11 

^12 

ll3 

12 

^22 

•^23. 

15 

■^3 



(3.4) 



47 


where A 22 ©’tc. are given hy the matrix equation, (3.4). 

The expression K[^u, v, wj becomes simpler for crystals of 
higher ssmiiretry and for particular reciprocal lattice points 
and particular directions of propagation of the thermal waves. 
Simple expressions for both theories of elasticity (using matrix 
Eq. 1.10 for old theory and Eq. 1.11 for the new theory) are 
derived in sections 3.1.1 and 3.1.2. 

The term for the volume of the crystaA (6?) is not occur- 
ing in the expression (Eq. 3.1) used for the diffractometer 
method for the reason that the direct beam intensity, is 
independently obtained in this method, unlike the photographic 
method (Eq. 3.2)as explained in detail in Sec. 4.5.5. 

3.1.1 Classical theory: 

Potassium Hiobate belongs to orthorhombic Bmm2 (Katz 
and Megaw, 1967) crystal class. The elastic constants in this 
system are given by the following matrix. (Krishnan, 1958, 
Bhagavantam 1966, Prasad and Wooster, 1956. Wooster, 1962). 
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Now from the matrix equation (1.10) we can get expres- 
sions for A. . for a general case as follows (Ramachandran and 

-L J 

T/Yooster, 1951 )s 


All 

11 

-h 

^ =66 

4* 

W'^Cr-c 

55 

+ OvwC^g 


2wuC.j3 

4* 

2uvC^g 

A 22 = 

"^'^66 

+ 


4- 

-4- 

0 

(M 

+ 2wwC2^ 

4- 

2 wuC^g 

4* 

2uvCgg 

11 

U^C[-r- 

55 

+ 

""=44 

4- 

w2c,3 

+ 2vwC^3 

4- 

2vnx0-z !~ 

35 

+ 


A 23 - 


4* 

’'^°24 

4* 

wS4 

+ vw(C23+0^^ 

)+vra(C^ 3 +C 3 g 

)+uv(C^g+C 


^12 = ^^^16 + '^^^26 + "^^*=*45 vw(Cg^+C25)+wu(C3g+G^4)+uv(C^2+%6^ 


...(3.6) 


With the help of elastic matrix Sq. (3.5) for D2» 

^2h crystal class the expressions A^^, A 22 etc. are simplified 
to the following equations; 


All 

= '1^0^ 1 + 

+ 


A22 

~ ’^^'^66 ■*■ ''^^22 

4 - 

*^°44 

A33 

= 

4 * 


C\J 

T- 

= (C^2 + Ogg) UY 



Ai3 

~ ^*^13 ^ 55 ^ 



^23 

= (C23 + C^^) TW 




With the help of these A,-- 5 's and the Eqs. (3.3) and (3.4) the 

J 

expressions for different K[^ u, for some directions 
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of propagation of thermal waves have been derived and are given 
in Tables 3.1, 3.2, 3.3. These relations have been used for 
evaluation of the elastic constants of Potassium ITiobaf e. 


3.1.2 hev/ theory of elasticity: 

The elastic properties of crystals belonging to D 20 

Dgjj) crystal class are given by the follov/ing matrix (following 
the notation of Laval, 1951). 


C^i 0^3 

"22 ^23 


0 0 


C 44 


77 0 

f 


^55 °58 


0 0 


0 0 0 0 0 0 

0 0 0 0 0 0 

« * ^ ~ ^ 


0 0 
0 0 


0 0 

0 0 

0 0 

t t 

^66 *^69 


(3.8) 


Now from the matrix equation (I.II) we can get the expres- 
sions for A . . for a most general case as follofire; 

J- J 


All ^ 


^ '"99 
55 


2rt ^ 

' ^66 

^22 

9 ? 

rOrjj 


^ ^88 
dcL 

dOjj 


2vw0g5 

"f- 

2uwC!|q 

+ 

2uvGi'g 

2vwC2^ 

+ 

2uw0jg 

+ 

2uvC2g 

2 'VV 70 ^rj 

-f 

2UWC33 

■f 

2 uvo!j^ 


contd 
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^23 " "^^^59 +"^^^34 +vw(C23+C4^)+uw(C^g+G^5)+uT(C^^+G^5) 


23 59 27 34 ^ ^'"'^^23 47 "^ ^ 39 45 ^ ^ 79 25 ^^ 

A ^2 — G^g + v Ggg +V 7 0^0 +vw( G 20 +G^g ) +uw( G^ ^+G 0 g ) +' u.v(G^ 2 ‘^Ogg) 
With the help of elastic matrix (3.8) for the particular case, 
(G 2 ^s I) 2 > ^ 2 h^ crystal classes, these relations are simplified 
to , 

^11 = ^1^011+ + "'^c'g 


122 

= 

+ v2c’ 2 


A33 

= uSo;^ 

2^« 

-h Y Orjrj 

+ w^G^^ 

A23 

= vw(G 23 

+ 0 ^^) 


•^13 

= uw( 0^0 

+ 0^^) 


■^12 

= uv(G^2 

+ Ogg) 



(3.9) 


Prom these relations (3.3) .and (3.4) tpe expressions 
for different kQ u, V jIaQ ^^^ 2 ^ for some directions of propagation 
of thermal waves have been derived and are given in Tables 3.1? 
3.2, 3.3 (Prasad and Wooster 1955). These relations have been 
used in interpreting the results of the present study in terms 
of the new theory of elasticity. 
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3.2 Evaluation of Elastic Constants 


3 . 2.1 Classical theory: 

It can be seen from Tables 3 . 1 ? 3 . 2 , 3.3 that constants 
Cii, C22. C^^, C55, Cgg can be determined independently 

from the observations along [j 0 ^ "'^001’ 

11,0 0 J okO L 0 1 0 3 001’ Co 0 3h00 E 0 OHoOl’ ^ ^hOO 

or j 100 Iq^^q respectively. The constants C^^g? ^235 0^3 cannot be 
determined independently, but can be determined by substituting 
values of independent constants 0^^, ^22’ ^33 

pressions for j~+ » + "Jp ’ ® J hOO ’ E+ '-Jz ’ ^ J 001 

Jo, + » + 'YaELokO respectively. Following the above 

method absolute values of all the nine constants of Potassium 
Niobate crystal (Crystal class C2^ or mm 2 ) have been determined 
from the quantitative study of its diffuse reflections by photo- 
graphic and diffractometer methods. 


3.2.2 hew theory: 

The new (Laval’s) 
elastic constants for ,an 
But for all the dynamical 


theory of elasticity requires fifteen 
orthorhombic crystal class £mm 2 j or C 
methods six of these would combine 


2v* 


into pairs and thus we are left with only twelve d3mamical cons- 
tants' to be determined viz. *^22’ ^33’ ^44’ ^55’ ^66’ ^77’ 

Ogg, Cgg,(ch+069), oonstan-ts 

^33’ ^44» ^66 classical (old) theory can be easily 

identified with C^^, C^3, C^^, o’^ and Cgg of the new theory? 


I 

i 
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t f r 

the additional constants of the new theory ^335 ^95 

determined independently from observations along ITo 1 03 qq]_? 
pO Ol'lhoo’ respectively. The constants 

which occur in pairs viz. (^25 *^47^’ ^^12 ^58^ 

Can be determined by substituting the above independent constants 

^ "1 -H1.L 

etc. into the observations along po, + ‘-^2 ’ -“ygJokO’ 

17+ PDhOO T? ’ - "T2^001 ^®sp actively. 

These constants cannot be separated out by the dynamical methods 

(Wooster 1961 ). 

In this way it is possible in principle to determine the 
absolute values of all the twelve dynamical- elastic constants of 
the Potassium Niobate crystal belonging to mm2 orthorhombic 
class, from the study of diffuse reflections only. The above 
derivations have been employed in interpreting and discussing the 
results of the present investiga.tions on Potassium Niobate in 
terms of the new (Laval’s) theory of elasticity in Sec. 5.1.2. 
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Direction cosines, 
u,v,w of the direc- 
tion of propagation 
of the thermal wave 
through node (h,k,l) 


TJiBLE 3 . 1 

a*h* Plane 

Expressions for 

K [ 

According to 
Classical Theory 


According to 
■''cw 

Theory 


L '' ° ^-IhOO 

C ° H OkO 

LO 103^00 
Ei jh’ -i °lh00 
Ui: 1/^’ i: i/T ^Hoko 




^(^22 ^66) 
2(C^1 + Ggg) 


C 


1 _ 

! 

11 

J- 

'99 


G 


66 


2(122 


99- 


2(1^^ + 


66 ' 


= LtfCag + °66^ + °22°66 - 2h2°66 " 

“ C + °66) °99^ " ^°-l2 * °69'>^ C 
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TABLE 3.2 
■b*c* Plaxie 


Direction cosines, 
u,v,v/ of the direc- 
tion of propagation 
of the thermal wave 
through node (hjk,!) 

Expressions for 

According to 
Classical Theory 

According to 

Nev/ 

Theory 

Co-ioJoko 

1 

^22 

1 

LOO^HokO 

1 

^44 

1 

"^44 

C ^ ^ 3 001 

1 

°44 

1 

C'rjrj 

L J: 71 ’ i irklioko 

2(044 °33^ 

2(c^^ + 

^2 

p ’ 

^2 

t - fk’ - 7^3 001 

2(*^22 ^44^ 

2(G22 + 0^4) 


f 


^2 ” rc !22*^^44 ^44^33 ~ ^^23^44 “ *^23^ 

^2 ~ r _ ^*^22 *^44^ ^^77 '*' *^33^ ~ ^*^23 ^47^ ^ 
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TABLE 3.3 
a*c* Plane 


Direction cosines, 
TijVjW of the direc- 
tion of propagation 
of the thermal wave 
through node (h,k,l) 

Expressions for 

K [ 

According to 
Classical Theory 

According to 
Uew 

Theory 

L ^ °3o01 

1 

C55 

1 

C55 

L lOO^hOO 

1 

"11 

1 

too^looi 

1 

O 33 

Olr, 

33 

L 0 ^ O hOO 

1 

^5 

1 

^88 

r + Y^» i T^looi 

2 (C^^ + 

+ °88) 

P3 


Cl 0’ ^T^lhOO 

2 ( 0^3 + 

2 ( 0 ^^ 4- 


^3 ~ C '^ 1 ' 1^^33 *^ 55 ^ ^ 3^55 ■” ^^ 13%5 ” ^13 ^ 

^3 “ L *^88^ ^%5 ^33^ " ^^58 "** ^13^ ^ 
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4. EXPEEIWTAl PROGEEUKE 


4. 1 Specimen Preparation 

Small crystals of Potassiiim Riobate ^parallelepiped in 
shape, were produced from binary melts (of K 2 C 0 ^ and Ebg^^) "^sing 
excess of £ 2 ^^^ as a flux (Wood 1951). Crystals, light blue, 
colourless and slight yellowish in colour were obtained with well 
developed faces and with a preferred cleavage perpendicular to 
Qo 10 ]] axis, larger crystals were grown by the technique des- 
cribed by Mller (1958). Crystals so grown from the above methods 
however have very large electrical conductivity near the transi- 
tion temperature (708°K) which makes the dielectric measurements 
difficult. Nearly single domain crystals of suitable dimensions 
for the diffuse scattering study were selected from them by a 
careful examination under a polarizing microscope, 

The volume of the crystal was determined either from the 
dimensions measured by a travelling microscope or from the weight 
and the density of the crystal. The second method involves less 
error, if the crystal weighs more than 300 jU, g (Burbank 196 5). 

Few crystals were finely ground to a cylindrical shape (diameter 
around 0.5 mm), one each for photographic and diffractometer 
methods. Crystals were annealed at about 800°C and were examined 
by laue photographs and under a polarizing microscope for any 
strains etc. 
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4.2 Comparison of Photographic and Coimter Methods 
(Wooster, 1962, 1964) 

There are essentially two types of experimental arrange- 
ments, one based on photographic and the other on ionization • 
(G-eiger, proportional or scintillation) counter, recording of the 
reflected X-rays. 

The apparatus required for photographic work is (i) a 
Eotation-Oscillation and laue Camera, with radius of at least 
5.7 cms. required for the accurate measurement of angular rela- 
tionships which are important in determining the direction in 
which the diffuse X-ray beam has travelled, (ii) A microdensito- 
meter for the measurement of the intensity of diffuse spots. 

The apparatus required for the ionization method is 
(i) a diffractometer, (ii) single crystal orienter, which is the , 
same as used in the single crystal - structure studies. 

The following differences between the photographic and 
ionization methods may be noted: 

(a) Even the weakest reflection can be recorded on the photo- 
graph by giving long exposures, unlike the diffractometer where 
the results on any single observation has to be recorded in few 
minutes; 

(b) In the photographic method a large volume of the recipro- 
cal space is recorded simultaneously, enabling a general survey 
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wliicTa is important in selecting proper planes for detailed stu- 
dies; the diffractometerjon the other hand, covers very small 
volume for a single observation. 

(c) Stability of the incident X-ray beam is very good in the 
photographic method because of long esposure and simultaneous 
recording of diffuse and direct beam intensity; while in the 
diffractometer the intensity of direct and diffuse beams are 
measured separately. This requires greater stabilization of the 
X-ray source. 

(d) Angular resolution is roughly 100 times higher in the 
photographic method than in the diffractometer method. The angu- 
lar resolution in the diffractometer can be improved at the 
expense of the number of counts. 

(e) The accuracy of intensity measurement of X-ray reflection 
by photographic means is hardly ever better than 5% (Wooster, 
1964) and to achieve even this accuracy requires very careful 
work. The diffractometer, on the other hand, can be arranged to 
give an accuracy of 1% with the help of very good stabilization, 
in most cases. Besides this^ a diffractometer is very convenient 
for measurements along given directions accurately passing through 
reciprocal lattice points. This property is very important in 
the determination of elastic constants. 

The photographic and ionization methods are each parti- 
cularly valuable in certain types of work. Each of them have 
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certain advantages and disadvantages. So it is difficult to 
decide which is hetterjin fact both the methods are in many 
ways complimentary. But in the elastic constant determination''-the 
ion isatioh' method is preferred, since it depends for its accuracy 
primarily upon the number of counts per unit time which can be 
measured in any given direction. Owing to random fluctuations, 
the number of counts which must be measured in order to attain 
an accuracy of 1% is 10,000. With the modem dif£cactometers 
10,000 counts can be registered in a reasonable time and hence 
no difficulty in attaining an accuracy of 1% (Arndt and Willis, 
1966). 

4.3 Photographic Method 

4.3.1 Preliminary adjustments; 

In the present investigation, Laue photographs were taken 
on a universal Rotation-Oscillation and laue Camera of 5.7 cm 
radius, designed and fabricated in our laboratory. The complete 
assembly drawings are given in Pig. 4.1 and 4.2 and the experi- 
mental arrangement is shown in Pig. 4.3. The cylindrical camera 
ha6 a collimator of 7 cm length. A conical rod resting on three 
balls fitted in a steel cup, passing through the axis of rota- 
tion of the camera carried a goniometer head which had usual 
arrangements for making axis of rotation of the crystal parallel 
and coincident with the axis of rotation of the camera. The 
rotation of camera could be read on the vernier attached to a 




PLAN 
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graduated disc, with an accuracy of 3*. A screw ca^tch carrying 
a vernier arrangement could fasten the rod carrying the gonio- 
meter head at any position, so as to keep the orientation of 
the crystal with respect to the incident beam fixed at any de- 
sired value. 

A crystal of suitable size was mounted on the top of a 
glass fiber (of diameter much smaller than that af the crjj'stal, to 
avoid any absorption of intensity due to this), with the help of 
an adhesive (Araldite or Quick Pix)? the whole thing was mounted 
rigidly on the top of the goniometer head with the help of a 
brass pin and sealing wax. Some aluminium powder was made into 
a thick paste in collodion and a small amount of it was put in 
the centre of the top of the crystal to give powder diffraction 
lines of AL superimposed on the laue pattern of the crystal. 

‘ These lines were used for (i) determination of the correct orien- 
tation of a particular plane with respect to incident beam, 

(ii) the measurement of camera radius and (iii) the calibration 
of angle of diffraction, jZi ? in laue pictures taken with mono- 
chromatic X-ray beam. 

After the preliminary adjustment was done optically, 
it was checked by taking rotation pictures (Appendix I, Plate 1 ), 
The faces being known, different directions in the crystal could 
be identified easily with the help of diffractometer technique, 
and so the crystal could now be rotated through a desired angle 
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orienting it in such a way as to record diffuse reflection from a 
particular plane (hkl). Finer rotations could he given with the 
help of a guide and a cam after fastening the rod with the catch. 

4.3.2 Orientation (laue) Photographs: 

Since we are interested in small wave vectors only, the 
crystal was oriented such that the incident heam makes an angle 
little different from the position of the Bra.gg reflection of the 
corresponding plane. After the crystal plane has been oriented 
the camera is adjusted in the X-ray beam, the tube being operated 
at 40 kV and 20 ml (Mo radiation). Because of the presence of 
white radiation, the whole Laue pattern will be recorded and 
because of the characteristic radiation, the diffuse reflection 
and powder diffraction pattern of A1 will be recorded simultane- 
ously on the film. Typical photographs are given in Appendix I 
(Plate 2). The indexing of the plane responsible for particulam 
diffuse reflection was ascertained with the help of a reciprocal 
lattice net and circle of reflection (orientation of the crystal 
being approximately known) and verified by diffractometer method. 

Since angles of diffraction, 0,for Al planes are known 
accurately, the distance on the film can be calibrated in terms 
of angle 0. Thus by measuring the distance between the A1(111) 
line and a particular laue spot the angle between the incident 
X— ray beam and a particular axis of the cr 3 rstal can be computed. 
The measurement of distances was done with the help of a 
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comparator and orientations thus determined are believed to be 

t 

ac curat e upt o 1 . 

4.5.3 Monochromatic Photographs: 

Por the determination of u, which are ultimately 

related to elastic constants it is necessary to determine the ratic 
of the diffusely scattered X-ray intensity by a particular wave 
vector to the incident intensity Iq. This ratio will be meaning- 
ful only when monochromatic radiation is used, (otherwise the , 

diffuse pattern due to one wavelength will be disturbed by the 
others) . 

In the present investigation MoKa (0.71071) radiation was 
used, operated at 35 KTP and 25 mA, 

The experimental arrangement used for taking laue pictures 
with monochromatic X-rays is shown in Pig, 4.3. Por monochromar- 
tization a plane crystal monochromator using a lithium fluoride 
crystal as a reflector was used, G. E. GA-»7 sealed Mb tube was 
used as a source of X-rays. The X-ray beam was first collimated 
through the rectangular slits S^, Sg and then allowed to fall on 
the lithium fluoride crystal mounted on a plate B, provided with ; 
necessary arrangement to have motion in two mutually perpendicu- 
lar directions. At first the crystal face was made parallel to ; 

the slits S 2 and then by rotating the crystal with the help of I 

a screw P^ it was brought to Bragg reflection position for (2,0,0)| 
plane (Bragg angle for MoKa -is 9 = 10.15°). The reflected beam 
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was seen on the fluorescent screen and the IIP crystal was so 
oriented that the reflected beam was uniform as checked photo- 
graphically. Spurious radiations at the edges of the reflected 
beam were eliminated by using straight knife edges and S^.The 
horizontal and vertical divergences of the monochromatized beam 
were also measured by exposing two films 11.5 cm apart. The 
divergences were calculated in the usual way by the measurement 
of the areas blackened in the two films. The horizontal diver- 
gence was found to be 8’ and the vertical 1°12'.' The latter 
assumed a much smaller value (22’) when collimated through slits 
of the camera and both are further reduced effectively due to 
small dimensions of the crystal. 

After the crystal has been oriented for a particular 
diffuse spot and an orientation picture is taken, the camera is 
collimated in the monochromatic beam for recording the diffuse 
reflection for intensity measurements. The ratio of the intensity 
of diffuse reflection to the intensity of the incident beam is of 
the order of 10“ . Since the intensity measurements by the photo- 
graphic method are reliable only if the blackening produced on 
the film lies (fall^) in the straight portion of the intensity 
versus blackening curve, it is necessary to reduce the intensity 
or the time of exposure for the direct beam such that the blacken- 
ing that it produces on the film is almost the same as that pro- 
duced by the diffuse reflection (both impressions being simulta- 
neously developed). The procedure followed in the present 
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investigation is to interpose suitable number of zirconium foils 
(.0002" each) of known thickness in the path of the direct beam 
just before it strikes the film. The thickness of Zr absorber 
is adjusted such that intensity of direct beam left over after 
absorption is nearly equal to the intensity of diffuse reflection. 
Uow, the reduced direct beam and the diffused beam both are allowed 
to fall on the film simultaneously and the film is exposed to both 
of them for the same length of time. The greatest advantage of 
this method is that the fluctuations of X-ray output of the tube 
do not affect the relative intensities, 

Por interposing Zirconium foils a small pocket of thin 
paper was made in the centre of the black paper envelope in which 
the film for taking monochromatic picture was wrapped. In this 
small pocket, suitable, number of Zr foils could be placed. As 
the intensity of diffuse reflection varies with the degree of 
misorientation and the structure factor of the plane concerned, 
the thickness of the absorber was varied accordingly, so that 
the intensity of reduced incident X-ray beam was nearly equal 
to that of the diffuse intensity. In the present study the number 
of foils were varied from 19 to 2l (i. e, the thickness varies 
from ,0038" to .0042"), It may be mentioned here that the absorp- 
tion coefficient of a substance is different for different wave 
lengths and in the monochromatized X-ray beam higher harmonics 
will be normally present. In order to avoid any error due to 
differential absorption of , V2, V3 etc. higher harmonics 
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were eliminated by running the Mo-target at 35 kV and 25 mA. 

This increased the exposure time considerably, which varied 
feom 35 hrs. to 50 hrs. in different cases. 

Each photograph thus taken would contain the diffuse ref- 
lectuion, the impression of the direct beam and the powder diffrac- 
tion lines of aluminium. To stop the back reflection lines (if 
any) of Zr from being recorded on the film, a thin sheet of lead 
was made in the form of a 2mm high, lidless and bottomless box 
which was hung from the edge of the cylindrical film holder or 
in the place of beam stopper so as to surround the Zr foils from 
all sides. (Typical monochromatic pictures are shown in Plate III,, 
Appendix I), 

4.3.4 Intensity Measurements: 

Following the method described by Sen and Chakraborty (1958) 
the intensity of the photographic records was measured using Moll 
recording microphotometer (Kippen Zonen), The film coiild be 
mounted on the sliding carriage of the photometer in the desired 
way. The scanning beam of light traverses the photograph when 
the carriage is moved. Since only the diffuse reflections in 
the equatorial layer line were considered in the present investi- 
gation, the photographs were mounted such that the scanning spot 
crossed the middle of the shadow of the crystal in the direct 
beam impression and the middle o'f the diffuse spot when the 
carriage was moved. 



70 


Tho dimensions of scanning spot were adjusted at 0,02 to 
,02 cm in height and 0.005 cm in breadth, the current for illumi- 
nating lamp was supplied by 6 volt accumulator batteries and was 
kept constant at 3.1 amp. The beam of light is focussed on the 
film, which has to be scanned. After falling on tho film., it is 
again converged on a vacuum thermopile. The current developed is 
recorded by a Moll sensitive galvanometer which is mounted on a 
vibrationless pillar. The time of response of the combination is 
of the order of 0.25 sec. The mirror attached to the galvanometer 
fibre reflects a fine beam of light which defLocts according to 
the density of blackening i.e. according to the current generated 
at the junctions of the vacuum thermopile and after being reflected 
by a total internal reflection prism is recorded on a strip of 
photographic paper. The paper is wrapped round a drum, the rota- 
tion of which is coupled with film translation by means of a 
variable gear. The movement of the reflected beam is parallel ’ 
to the axis of rotation of the drum. The gear ratio used in the 
present investigations was 1:4. The distance between the two peaks ' 
corresponding to A1(111) lines in photometric record was used in ; 
converting the distances (abscissa) of the photometric curve into i 
corresponding angular coordinates of the photograph. Tho intensitiej 

were obtained by comparing the ordinate of the photometric curve i 

1 

with the standard wedge following the method of Robinson (l933).Th! 
calibration wedge was recorded for each picture separately on the pj 
tographic film cut out of the same sheet as used for recording the | 
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mono cTiro mat ic picture and both were developed together. The time 
of exposure is linearly proportional to the distance along the 
length of the wedge and hence the intensity corresponding to any 
ordinate of the photometric curve is proportional to the dis- 
tance from the zero intensity point. The photometric curve for 
the wedge was recorded under the same photometric conditions as 
the diffuse spot and the direct beam impression. Typical photo- 
metric curves for the calibration wedge (Appendix III^C) and 
also for the diffuse spot and direct beam in 5 )ression have been 
reproduced in Appendix The intensities of diffuse 

reflections and the incident beam were thus measured from a photo- 
graph on the same arbitrary scale. 


4.3.5 Evaluation of relative intensity- : 

The direct beam impression will have in its centre 
the shadow of the crystal due to the absorption by the crystal, 
hence its micro photometer record will have a dip in the centre 
corresponding to the shadow of the crystal (Appendix III, a). 


The intensity corresponding to the ordinate for the dip 
of the photometric curve of the direct beam impression (the 
reason for taking the dip has been discussed in Appendix III) 
may be denoted by I^ i.e. reduced intensity of the incidoQt 
beam at the film. 


'0 


. M-t 


A 


Eow 
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is the intensity of the direct hGam at the film, 

^ the linear absorption coefficient of Zirconium and t^ the 
thickness of the zirconium foils used for reducing the intensity. 

t 

Prom we can compute the intensity of the direct beam at 

the crystal, let the area over which the direct beam spreads at 

the crystal and on the film be S- and S-m respectively. Then it 

S ^ s 

is evident that i = i ^ . This factor ^ can be evaluated 

° ° ^c ^c 

from the knowledge of the divergence of the beam and the size 
of the slit. In the present case it was 1.28. Next comes the 
consideration of the value of linear absorption coefficient of 
the zirconium foil for MoKa, which w<ls used to reduce the inten- 
sity of the direct beam. In order to avoid any error cansed due 
to variation of the purity of the sample of zirconium foils the 
value of mass absorption coefficient was determined by the 

C st«- o-ppeTid/x .zzr) 

diffractometer method as described by Chakrraborty (l958)i’ It 
was in full agreement with the one given in the International 
tables for X-ray crystallography (1963) which shows that the sam- 
ple is highly pure. 

The intensity corresponding to the ordinates of the photo- 
metric curve of diffuse reflection is found out from the photo- 
metric curve of the calibration wedge. The distance betv\feen 11 
(ill) peaks in the microphotometer record was calibrated in 
terms of the angle of diffraction. Thus the abscissa of photo- 
metric curve for diffuse reflection gives 0 and the corresponding 
ordinates give the observed diffuse intensity l^ (uncorrected 
for any factor) . plotted against 0 for each case and 
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these Curves were used for determining the value of I(j/lo corres- 
ponding to a certain 0 , which is rela,ted to a particular wave 
vector and was found out by calculation. Some of these curves 
have been reproduced, shown in Pig.l- 14 , 

4.4 Diffractometer Method 

Since the diffractometer with single crystal orienter 
(Goniostat) permits a continuous survey of a whole hemisphere 
of a reciprocal lattice with smooth and finely controllable 
motions, it is ideally suited for holding and manipulating crys- 
tal specimens for a detailed study of their diffuse scattering. 
Without changing the physical arrangement and alignment of the 
Goniostat and SPG spectrogoniometer (Pumas 1957) as it is used 
in the single crystal structure studies, remarkably good results 
can be obtained through the use of balanced filters (for mono- 
chromotization) for the X-ray beam. The absolute value of the 
direct beam intensity (1^) is obtained from Compton scattering 
of diamond as intermediate standard (Appendix II). 

4.4.1 Method of Monochromatization: 

The system of balanced filters (Ghipman and Warren (l950)j 
Tanaka et al(l959),, Guinier (1945), Peiser et al ( 1955) ,Umanskii 
(i 960 ), Burbank (1965), Arndt and Willis (1966) and Ross (l928)) 
is an alternative to the monochromatizing crystal. The filters 
are precisely matched so that they transmit the same proportions 
of the radiation on either side of the small spectral range which 
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contains the Ka line of the given X-ray target. The advantage 
of this method of monochromatization is that the intensity of 
the monochromatic beam is much higher than it is with crystal 
reflected radiation. Unfortunately, there is small portion 
of the white radiation which is transmitted more by one filter 
than by the other, but this usually only involves a very small 
correction to the results, 

4.4. 1.1 -Balanced Bilters: 

The operating principle of balanced filters (Ross 1928) 
is shown in Big. 4.4. Two filters are constructed, generally 
of two elements of atomic numbers Z and Z- 1 , whose absorption 
matches at all wave lengths other than over the narrow band 
between the two absorption edges X 2 2-1* mea- 

surements are made, one with each filter, the difference between 
the two is due solely to radiation within this narrow ‘pass- 
band’. i.e. monochromatic radiation. 

However in the present case we have used Zr and high purity 
'^ 2^3 ^ pair of balanced filters for Molybdenum target. 

The latter filter is used in the form a thin pellet of the 
oxide of the metal because pure yttrium metal in the form of 
a thin foil was not available. The balahbh obtained for (004) 
reflection is shown in Big. 4 , 5 . 

The balance of the filter pair must be tested at every 
stage by plotting transmitted intensity as a function of 
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the wavelength for both filters, The wavelength scan can he 
performed hy moimting an analysing crystal on the diffracto- 
meter and hy varying 26 and 0 (SPG- goniometer motions) in the 
ratio 2:1 as can he seen in the Pig. 4.5. 

4.5 Correction Pactors 

The crystal used on the diffractometer is fully immersed 
in the X-ray heam hy using a suitable beam collimator; hence 
many of the corrections involved can he applied in a similar 
way as in the photographic method. 

4.5.1 General Scattering: 

In the diffraction process hy the diffractometer and the 
photographic technique a general scattering underlies Bragg 
and diffuse scattering. The general scattering includes, Compton 
scattering, flourescent, radiation and scattering due to optical 
waves. The intensity of the general scattering varies to a 
limited extent over the region of the reciprocal space, covering 
the observations about one relp. 

The correction to he applied to the observed diffuse 
intensity due to general scattering is the same in both the 
techniques. 

In the diffractometer method such variation in the 
general scattering can be eliminated by talcing the mean value 
of the intensities of X-rays scattered from reciprocal volumes 
at equal distances on opposite sides of the relp. The diffuse 
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intensity due to first-order scattering from elastic waves, 
is found by plotting observed diffuse intensity, (corrected 
for skew, polarization, absorption and divergence) against 
(i.e.}\^). The points so plotted lie, within the accuracy of 
the measurement, on a straight line, which cuts the axis of 
at a point corresponding to the const aht contribution of the 
Compton, fluorescent and optical scattering. This quantity 
is subtracted from the observed intensity to obtain the inten- 
sity due to the first and second order elastic scattering. 

In the case of the photographic technique, intensities 
due to diffuse reflections are superimposed over the background 
intensity due to Compton scattering, air scattering etc. By sub- 
tracting the extrapolated backgroimd intensity from the total, 
the corrected intensity for the diffuse reflection only has been 
obtained. The background intensity has been shown as a dotted 
line in the 0 curves in Big, I-I 4 (Appendix I). Thus 

the values of given in the table A-I 6 (Appendix I) are 

corrected for general scattering. 

4.5.2 Divergence correction: 

An ideal collimator gives rise to a perfectly parallel 
beam of X-rays. However, the finite size of the slits between 
the crystal and the X-ray tube and also between the crystal and 
the detector permits X-rays having certain divergence to enter 
the detector and be counted together. The experimental 
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observa-tion does not, therefore, correspond only to the intensity 
of diffuse scattering from a single point in reciprocal space 
hut also covers a certain volume of reciprocal space surrounding 
the point, The greater the divergence of the X-ray heam the 
greater the volume of reciprocal space over which the integra- 
tion is effected. 

4. 5. 2.1 Photographic method: 

Phe divergence of the incident heam in this case is very 
low (horizontal 8' and vertical 22*) because of the very narrow 
slit employed and the short distance between the crystal and 
the collimator. ( aj 2 cm). The effective divergence would be 
further reduced due to the small size of the crystal (0.996 x 
10”'^ cc). The divergence correction was determined in few cases 
and was found to be negligible, 

4. 5.2. 2 Diffractometer method; (Wooster 1962) 

The size of the collimator slits and the crystal used in 
the diffractometer method is comparatively larger than in the 
photographic method. As a result the divergence correction be- 
comes appreciable. However in the present study unusually 
small crystals were used. The total divergence may be divided 
into three components denoted respectively by i, 9 andX* 
i-divergence is the maximum angle between the rays falling on 
the crystal when they are projected on a horizontal plane. The 
6 divergence is the maximum angle between the reflected rays, 
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also projected on a horizontal plane, which can enter the 
detector at any given setting. The 3C 

maximum angle between the reflected rays, projected on to a 
vertical plane passing through the crystal and the detector, 
which can be recorded at the given setting. The magnitude of 
divergence will vary from one experiment/ to another, depend- 
ing on the size of the crj^^stals and the slits; in the present 
experiment i is about 30', 0 about 1° and ^ about 5°. 

In order to apply the i-divergence correction, the de- 
tector and X-positions are kept fixed. The intensity of the 
Bragg reflection is measured as a function of the angle i in 
the vicinity of the Bragg peak position (Big. 4.6a). The angle 
i is read on the 0 circle of the crystal orienter. Then inten- 
sity data (I)(i)) are collected at a fixed positions of i (such 
as A in Big. 4.6b) which are atleast l'^ away from the i corres- 
ponding to the Bragg peak. The angle i is varied in small 
intervals (6 ir'>0,02°) around point A, to obtain the curve AB 
in Big. 4.6b. The diffuse intensity varies slowly with i, where- 
as the Bragg reflection varies rapidly with i in the neighbour- 
hood of the Bragg setting (Big. 4.6a,b ), Each point on the 
curve AB can be regarded as being produced by many planes sligfetl-y 
slightly inclined to the reciprocal lattice planes which give 
the Bragg reflection. At a point A, Big.4.6b,one lattice plane 
is in the Bragg setting and its contribution to the intensity 
is represented by the small copy CD of the Bragg intensity 
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curve drawn at the same angle i as that of point A. This parti- 
cular plane contributes the maximum possible diffuse intensity 
at point A, Other planes nearly parallel to the A set, contri- 
bute an amount corresponding to the departure of the angle i 
from that value which corresponds to their maximum contributions. 
Thus at any point A, the value of D(i) is the sum of contributions, 
each one of which is the product of the maximum intensity of dif- 
fuse scattering from a lattice plane, and an ordinate of the curve 
GD.Ie''ti’.B denoti^by D(i) the true intensity curve of the diffuse 
reflection due to the disturbed lattice without any modification 
due to the shape of the B(i) curve. The angular breadth of the 
B(i) curve is denoted by , and the angular separation of any 
point on the B (i) curve from the central maximum is called 5 
(Big. 4.6b). The particular value of i at the point A on 
curve D(i) is denoted by i^. The elastic waves corresponding to 
a Bragg setting (i^ +6) give a contribution at A eq.ual to 
D*(i^ + 6) B(6), where B(6) is the ordinate of the B(i) curve at 
a distance 6 from the maximum. The total effect is thus the 
integration of all these contributions from the separate elastic 
wave trains, namely, 

+ A/2 , 

DCi^) = / D (i^ + 6) B(6) d6 (4.1) 

- £^/2 

It is assumed here that the Bragg reflection curve is symmetri- 
cal about its centre. This assumption is usually justified. 
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If the Bragg reflection curve is not symmetric with res- 
pect to its centre then the error due to this assumptior is 
removed hy taicing the values of B(i) at + 6 and - 6 (Pig.4.6c). 
Brom an experimental curve such as Big. 4.6c, which gives I)(i) 
as a function of i, we may calculate the quantity &(6) for each 
value of 6 Between 0 and + ^/2 according- to the expression 

G(6) (I)(iQ+6) + D(i^-6)) - ^^*2) 

The expression G-(6) gives the correction factor at every value 
of 5 which must he applied to the experimental curve, lo take 
into account the shape of the B(i) curve, each value of G(6) is 
multiplied hy the corresponding B(6) and integrate the product 
over the range + <6/2. This is divided hy the integral of B(6) 

(to normalize) over the same range to obtain the final correction, 

+ A/2 + h/2 

Z(i ) = / G(6) B(6) d(6) // B(6) d6 (4.3) 

° - A/2 - A/2 

The evaluation of this correction is carried out graphically, The 

observed value of the diffuse intensity at any given point i^ is 

multiplied hy a factor (1 + Z(i^)) to obtain the correct -value 

of diffuse intensity, 

0 - divergence and X -'ii’^ergence corrections are carried 
out in a manner similar to i - divergence correction. In the 
former i andX kept in a fixed position and intensity is 
collected as a function of 9, while in the latter i and 8 are 
kept constant and intensity is studied as function of X. The 
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intensity data are analysed in an analogous manner to determine 
the correction factors. Finally i, 0 and X~di'^63:gence correc- 
tions are added to obtain the corrected diffuse intensity. 


4 . 5.3 Skew correction: 

The skewr correction arises because the glancing angles of 
the incidence and of reflection are not eq.ual. laval (1939-51) has 
shown that diffused intensity in such a case is given by 


» -I 

- ^ "2 


1 + 


Sin i 
Sin (0-i) 


(4.4) 


Where is the observed diffuse intensity i is the angle of 
incidence and 0 is the angle of diffraction. The factor 
1/2 |l+ Sin i/Sin(0-i)J is called skew correction. The observed 
intensity must be multiplied by the skew correction to obtain 
the intensity which would be observed if the angles of inci- 
dence and reflection were equal. It will be seen that when 
i > 0/2 the skew factor is greater than unity and when i < 0/2 
it is less than unity. When i > 0/2 the voltune of the crystal 
irradiated is less than when the Bragg condition is satisfied (i.e. 
when i = 0/2). For both these reasons is less than This 

correction is the same in the photographic and diffractometer 


methods and has been applied to each observation. 


4.5.4 Polarization factor: 

The polarization factor is given by, 

P = (A^ + Cos^ 0) 


(4.5) 
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where, A and B have the same meaning as in Eq..1.l3. If a mono- 
chromatizing crystal is not used as in the diffractometerjthen 
the radiation is unpolarized and, 

A = B = l/f2 

A^ = l/l+Cos^2e, B^ = Cos^2e/(l+Cos^2e) 

9 is the angle at which X-ray beam is reflected. On the other 
hand, since the crystal reflected monochromatic beam has been 
used for recording all the photographs from which intensity 
measurements have been made, the polarization factor for equa- 
torial reflections tahes the form, 

p = 1 + Cos^ 2a Gos^ 0 

where, 2a is the angle of diffraction of the monochromatic beam 
plane and 0 the angle of diffraction at which the observation 
is made. A general expression applicable to non- equate rial ref- 
lections also has been given by Azaroff (1955). 

4.5.5 Absorption correction: 

let be the initial intensity of an X-ray beam incident 
normally on the absorbing material (crystal) and I the intensity 
after the beam has travelled through a distance 3c in the absorber. 
If dl represents a further diminution of intensity of the beam 
on its travelling an additional small distance dx, then-dl/dx 
is the rate of decrease of intensity. If it be assumed that 
dl/dx is proportional to I, we may write, 
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dx 


Ail 


(4.6) 


where Jk. > the constant of proportional it jr, is called the linear 
absorption coefficient, which is a characteristic of the absorbe 
crystal and is a function of the wavelength of incident X-rays. 
On integration and applying initial conditions, we get, 

I = e (4.7) 

If the distances traversed insic’e the crystal by the inc 
dent and the diffusely reflected beams are equal, then the inten- 
sity of both of them will be reduced to the same extent i.e, by i 

jjl^ IT / 

factor e" and consequently the ratio I^/Iq remains unaltered, 
Such will be the case (for equatorial reflections) when the crys- 
tal is spherical or cylindrical. This is the reason why yU. does 
not occur in the expression for intensity formula ( Hq. . 5.2) ' tVi 
the photographic method. Some times it is not possible to obtain 
cylindrical or spherical crystals because grinding the crystal 
may lead to fragmentation or twinning; shaping is also difficult ! 
if there are strong cleavage planes. Therefore, for an unground . 
crystal, (if the dimensions of the faces boimding the crystal are 
measured) the distances traversed bylhs 'incident ^nd diffuse 
beams would not in general be equal and correction might be neces- 
sary in certain cases, for applying this correction (in the phot( 
graphic and diffractometer methods) the value of^kshould be know; 
Simple methods for the determination of /^by the photographic and 
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diffractometer techniq.ues are described in Appendix III. A 
graphical method due to Albrecht (1959) for treating absorption 
in a crystal of any size, shape and absorbing power is briefly 
described here. 

In an actual crystal the path length x is a variable 
depending ion the size and shape of the crystal and on the dif- 
fracting plane. The intensity of the diffracted beam is given 

1 = 1/ dv/v (4.8) 

V 

where I^ is the intensity which the diffracted beam would have 
possessed had it not been for absorption. This integral is 
difficult to evaluate but an approximate solution can be obtained 
graphically by dividing the crystal into small volume elements 
and summing the absorption effects of each,Eq.. (4.8) may then be 
rewritten as, 

^ _ lo f dT _ e-^1 

V ^ 

where n volume elements are considered. In Eig. 4.7 the absorp- 
tion by a volume element is given by l/l^ = 

t^ + tg is the path length (in the crystal) of beam diffracted from 
this volume element. If the incident beam .Sq is divided into 
a number of equal rays as shown in Eig. 4.8, where 0 is the Bragg 
angle for the plane hid considered here, then S will represent 
the diffracted beam, divided into rays of the same dimensions as 
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HCr. 4,8s Projection of the crystal | showing 
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Sq. If the pTOjection of a volume element is represented by the 
intersection of the two rays, it is possible to obtain a value 
for the summation of the Eq.. (4.9) by adding quantities of the 
type and dividing the sum by n, the number of ray inter- 

sections associated with the diagram. The number of intersections 
within the boundaries of the crystal is equal to the number of 
volume elements chosen. The path length for a given volume ele- 
ment can be measured with rulers and the quant itjr can be 

readily obtained from a table of exponentials. It is thus only 
necessary to draw a cross-section of actual crystal to the scale; 
draw in the incident and emergent rays in the direction appro- 
priate to tbe plane under consideration; measure the path lengths, 
and count the number of intersections. If a large number of cal- 
culations are to be performed the work may be carried out on 
computer. In the present case calculations of absorption factors 
are carried out on IBM 7044 Computer, 

The amount by which the intensity of the hkl reflexion is 
reduced by absorption is denoted by The reci- 

procal of is the absorption factor ^ repre- 

sents tile factor by which the observed intensity is multiplied 
to obtain the corrected intensity. 

4.5.6 Second order diffuse scattering correction: 

Main contribution to the diffuse intensity is due to 
first order diffuse reflection. Prom equations 1.16, 1.19 and 1.25 



the ratio of the intensity of first and second order diffuse 
reflection is given hy 
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(X~' 5 

™ = “^ kT X q K C '^’"'^’''''^Ilhkl ^ K[] u, v,w'3 (4.10) 

In the photographic method used for (4 0 0) reflection the theimal 
wave vectors are ver 3 '' small and hence the correction is neglected. 
However in the diffractometer method correction has keen applied 
to (6 0 0), (0 0 6), (0 4 0), (4 0 0) and (0 0 4) relps: as suggested 
hy Ramachandran and Wooster (1951), Prasad and Wooster (1955) 
and Amoros (1968). 

f | m -ifn O 

.*. K C Jhkl ~ ^ C ^’"^’''''^Ifhkl (Approximation) 

*** ~ ^ ^ ^r.^>v>^Il hkl 

or 

3 

0^ ==(!:[ ^ kT K[ 

The diffuse intensity due to first order only is obtained 
hy subtracting CT^ from the observed diffuse intensity (jHj (correc- 
ted for divergence, skew etc.). Third and higher order corrections 
are fairly small and hence neglected. Yery recentlj?- a least square 
fit method for the deteimination of the above correction is given 
by Lucas (1968, 1969). This method is particularly useful when 
the second order correction is large and very accurate values of 
the same are needed. However in the present investigation appre- 
ciable difference between the two methods of correction are not 


found. 
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4,5.7 White radiation and mosaic - structure corrections: 

Unlike the skew correction and the divergence corrections, 
which have to he made for every measurement, the corrections consi- 
dered here arise only for measurements made at points along certain 
rekhas. Thus, if the crystal is appreciably mosaic, then the Bragg 
reflexion would extend over a range of angles on either side of 
the peak, e.g, along the arc BQB* in the section shown in Big, 4. 9. 
Thus, enhanced values would he observed for the diffuse flux for 
points lying on rekhas close to the rekha AQil* . If this additional 
flux is not large compared with the true diffuse flux, then it can 
he eliminated by interpolation from measurements along neighbour- 
ing rekhas. This procedure of interpolation is called the correc- 
tion for mosaic structure. 

The white radiation correction is important for rekha G'qC, 
i.e. rekha lying along the rel- vector OQ and is similar to mosaic- : 
structure correction, since it is intended to correct for a slightr, 
enhanced diffuse reflexion obsenved for points along the rekha. The ; 
enhancement is due- to small amount of white radiation diffusely 
reflected by the crystal under study. The enhancement is generally 
less than -10 percent, and it can be corrected by interpolation 
from measurements at points such as C 2 G^ rnd c lying on the 
same circle as but on rekhas at + 10° and + 20° to the rel- 
vector at other angles as suitable. It is however important to 
make this correction for every measurement lying on this rekha* 
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The final intensity ( 1 ^) is obtained after applying all these 
corrections. 

4.5.8 Trigonometric evaluation of a Thermal wave vector: 

The wave vector q (QP^ QP^, OQ-^ etc. 'Pig. 4.10) for any 
particular direction of propagation of the thermal YJ'ave and the 
corresponding angle of scattering 29 can he calculated hy estab- 
lishing trigonometrical relations by considering the geometri- 
cal position of the circle of reflection in the reciprocal lattice 
net for a particular setting of the crystal. The position of the - 
circle is completely determined by the knowledge of the correct 
orientation of the crystal with respect to the incident beam. 

This can be readily determined by measuring the angle of diffrac- 
tion of some known reflection in the orientation pictures. The 
parameters used for making these calculations are given below: 


* 

a 

= 0.1755 

^-1 

b* 

= 0.2518 


c* 

= 0.1748 

^-1 

1 

= l.^-O 

^_1 


'K* 'H' 

where a , b , c are imit vectors in the reciprocal lattice. 

It is evident from the equation (3.1, 3.2) that graph between 
and for a particular reciprocal lattice point and direc- 

tion of propagation of thermal wave should be a straight line 
passing through the origin (after subtracting general scatter- 
ing correction). The slope of this straight line would give 
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-*2 

q. . These graphs we 2 ?e drawn for all the cases and were used 
for evaluating (or in the photographic method). 

Ill other factors in Eqs. 3.1 and 3.2 being known the value of 
K;[^ u, have been calculated. These graphs have been re- 

produced in the Appendix I. 

4.5.9 Equiscattering Contours around a node and K surfaces: 

A surface surrounding a reciprocal point over which 
Id/^o ^ constant value is called a iso-diffusion or equi- 

dif fusion contour. 

The procedure for drawing equiscattering contours in 
the photographic method is discussed first. 

In order to see the variation of diffused intensity with 
the wave length in different directions of propagation of ther- 
mal waves responsible for diffuse reflections, equiscattering 
contours around (hOO) reflection in a*c* and a*b* reciprocal 
lattice planes can be drawn in the following way, A curve 
I^/lo vs. 0 for a (hOO) reflection in a*c* and a*b* planes 
corresponding to different orientations is obtained. Erom this 
curve corresponding to a maximum misorientation from the Bragg 
reflection position a sufficiently high value of I^/Iq (of the 
order of 10”^), was chosen, because this value will be available 
in all the curves. From each of the I(^/Iq versus 0 curves the 
approximate angles (0’s since each curve will give two angles) 
which correspond to the above chosen value of I^/l^ were noted 
down and then the skew and polarization factors for these' two 
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0’s were calculated. The required corrections thus being appro- 
ximately known, the values of which will give the corrected 

■value equal to the chosen value of at 0^ and 02 were easily 

found out for each case. 

An equiscattering contour is drawn on a reciprocal net as 
follows (Fig. 4.10). A reciprocal net on a large scale (0.1 
= 1 cm) was drawn. Using the direction of incidence of X-rays 
as the diajneter, the circle of reflection v/as drawn with, ra- 
dius = 1/ X where X= 0.7107 i for MoKa, radialion. Then, a 
line corresponding to the vector of diffusion (i.e. line making 
the angle of diffraction 0 with the diameter of the circle of 
reflection e.g. OQ^ in Fig, 4.10) was drawn. Intersection of 
vector of diffusion and the circle of reflection gives the corres- 
ponding pole of diffusion. Thus poles of diffusion corresponding 
to a particular value of (or in photographic method) at 

different orientations for a particular diffuse reflection were 
plotted in a reciprocal lattice plane. These curves are- -known as 
equiscattering contours or iso-dif fusion contours for a parti- 
cular diffuse reflection. 

Uow turning to the diffractometer method, the procedure 
for drawing equiscattering contours is exactly identical as 
described above for the photographic method, except that in the 
diffractometer technique the graphs for a given reflection are 
drawn between (corrected diffuse intensity) and X^» instead 
I^/Iq vs. 0 as in the photographic method. This is so because 
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and are measured separately in tlie diffractometer method, 
while they are recorded simiiLtaneously in the camerd, method. 

0 corresponds to a definite value of , the wave length of 
the elastic wave. A graph of vs. X ^ is a straight line, 
drawn on hoth sides of a given relp. Erom these tvro graphs 
of 1 ( 3 ^ "^s. A I for a given relp, a value of corresponding to 
maximum misorientation from Bragg reflection position such that 
a sufficiently high value of is obtained. Dividing this by 
it gives the value of Now substituting this value of 

I^Ao for a given relp in the intensity formula for the diffrac- 
tometer ( Section 3. 1. 1),i.e. , 




( ■=“ ) = Constant value, 

■^o hkl 


One gets, 


= Const. K[^ u,v,w3 


hkl 


(4.11) 


A curve can be traced which is an equiscattering contours 
for the given relp. Such equiscattering contours are drawn for 
(40 0) and (6 0 0) nodes (See. 5.1.4). 

4,5.9. 1 K-surface for Orthorhombic EZNbO^: 

The variation of the diffuse scattering from poirrts round 
a given reciprocal lattice point can be conveniently expressed 
by means of a surface, first used by Jahn (1942), Prasad and 
Wooster (1956) which we shall describe as E-surface. This is 
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defined by -Eq. (1.19). 

2Wi^^ + 2NLA^^ + 2LM“^ ... 

Thus the value of Kj^ u, v,w^ is determined simply by 

elastic constants and the directions of (a) the wave normal 
corresponding to the source of the diffuse scattering, and 
(b) the rel vector. The expression for K [u, Vy in terms 
of the elastic constants and the direction cosines contained 
in the equation (1.19) is applicable to any system of symmetry. 
As the symmetry decreases more and more elastic constants be- 
come involved and the labour of computation increases. By 
inserting the particular values of elastic constants and giving 
u,v,w all possible values, the whole range of K- values may be 
covered for a given relp, K-surfaces are drawn for (0 0 6) and 
(6 0 0) relps projected on [["O 1 0^ and [[OOl^i 3:espectively 
(Sec. 5.1.4). 

4.5.10 Choice of the planes for the determination of the 
elastic constants ■^**)* 

' The diffuse flux to be observed should be as large as 
possible, since the intensity of diffuse reflection in the first 
order depends directly on and i.e, 
tions remaining the same. Therefore, it is best to choose for 
study those relps for which this quantity is largest. The con- 
venient range of suitable Bragg angles is 20° - 50°. The lower 




T^^hkl > 


other condi- 


1^ A 


11 


+ 






•35 
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limit arises from the fact that the skew correction is appre- 
ciably different from unity for angles of reflexion smaller 
than 20°, and slight missetting of the crystal is liable to 
produce large errors. The upper limit arises from two causes 
(i) The mechanical limit to the movement of the instrument 
which goes only upto 120° conveniently, so that 0 should he 
less than 50°, (ii) The formula used hy laval, based on the 
assumption that the amplitudes of atomic vibrations are small 
compared with the lattice spacing. The mean amplitude of 
thermal vibrations at room temperature is of the order of 0.2i 
for most crystals (Lonsdale 1942), we may therefore assume 
that'/x>1 i is a safe lower limit for the lattice spacing. 

The Bragg angle corresponding to 1 i with Cu Ka is ilii.50°, 
while with Mo Ka Qi. 25°. Considering now the wavelength of 
X-rays to be used, it is preferable to use as long a wave length 
as possible, for several reasons as explained in detail ' by 
Ramachandran and Wooster (1951), . But the absorption factor 
has also to be taken into account. It is very large for a long 
wavelength , thereby causing a large error in intensity. Hence 
the choice of proper wavelength .for a given case should be made. 
Bor the planes with larger X (i.e. higher Bragg angle), second 
order diffuse scattering will be larger (as discussed earlier). 
Therefore correction on this account will be necessary; P, the 
polarization factor is smaller at higher angles and diffuse 
intensity will be reduced due to this factor. Hence, planes 
which have high Bragg angles can be used only if the values of 
their structure factors are very high. 
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TABLE 4 . 1 


A. Planes studied with 

0 1 0 '[] axis 

vertical 

Indices 

d in S. 

Diffract ion angle 
2 6 in degrees 

?! 

4 0 0 

1.424 


28'. 88 

72.73 

0 0 4 

1.430 


28.76 

71.92 

0 0 6 

0.952 


43.78 

38.77 

B, Planes 

with 0 0] 

axis vertical 


0 4 0 

0.996 


41.94 

62.83 

C. Planes 

with [0 0 

0 

axis vertical 


6 0 0 

0,946 . 


43.96 

41.44 


S'tniC'fcure factor dafa are faiceii from Katz and Megaw (1967) 
and Dwivedi and Srivastava (1969) • 
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5. RESULTS ML DISCUSSION 


5.1 Results 

5.1.1 Elastic constants; Old tUeory 

The diffractometer study of orthorhomhic Potassium Niobate 
crystal at room temperature resulted in the following nine elastic 
constants, based on old (Voigt’s) theory (Table 5.1). Various 
corrections, described in Sec. 4.5 have been applied to the dif- 
fractometer data before obtaining the final results. The details 
of calculations of these elastic constants are given in Appendix I. 


Table 

5.1 


Elastic constant 


Value 

^11 

0.5(54) 

i 9 9 

X 10 dynes/ cm or 



10"^^ Uewton/Meter^ 

CM 

CM 

O 

0 . 5 ( 16 ) 

ft 

C 33 

0 . 6 ( 32 ) 

tf 

°44 

3.5(58) 

ft 

°55 

3.6(51) 

If 

°66 

2.0(69) 

!f 

C 12 

-1.0(58) 

If 

O 13 

-2.0(27) 

If 

CM 

0 

-1.9(14) 

tr 
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5.1.2 Elastic constants: New theory 

It is seen in Sec. 3.1.2 that fifteen elastic constants 
are necessarj^ for completely specifyin^i the elastic behaviour 
of Potassium Niobate crystal. But six of them do not occur 
independently but appear only in combinations 
( 02 ^ + *^^-13 ^58^ expressions for propagation of 

elastic waves in the crystals for the pa^rticular case; therefore, 
all the dynam.ical methods would only give these combinations and 
not the individual constants (Wooster 1961). One is thus left 
with only twelve dynamical constants. Namely, 0^2? 

Crjrj,0QQ, Ggg, ( C ^2 *^ 69 ^’ ^*^23 ^ *^ 47 ^ ^^13 **” *^ 58 ^* 


Out of these 0^^ to can be identified with the constants C 


"66 


11 


to Ggg of the classical theory (fables 3.1 to 3.3). So the values 
of these constants are unaffected due to the implications of the 
new theory as applied to the present case. 

further, it is possible to distinguish between the follow- 
ing pairs of elastic constants in the new theory: (G^^, Gggj Ggg, 

G 99 5 G^^, Gyy? *^ 69 ’ *^66’ *"58’ *"55 ’ *"47^ *"44^* exan^ile, 

according to old classical theory the reldaas K[) ''GGHqqx 

G iHboO give the same elastic constant, namely G^^ 

but according to new theory, they give rise to two different 
elastic constants i.e. K[^ 100 []qq 2 gives G^^ while KQ00 1 J|^qq 
gives Ggg. Now it is known that is identified with 
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0^^ and Cgg laave been determined independently from indi- 
vidual valuer of lijj i 0 0 J and ICj^O 0 ^qq respectively. Details 
of observations and calculations are given in the Appendix I in 
the tables for (O 0 6 ) and (4 0 0) planes; whence we get, 

0^^ = 3.6(5) X lo”'^ dynes/cm^ 

C 38 “ 3 . 4 ( 6 ) X lo”^^ dynes/cm^ 

The difference between the two is about 6*^^. Since the two values 
depend on the structure factors of the two different planes and 
because of the limitation of the accuracy attainable by the method 
(Sec. 5.4) used in the present investigations, this difference is 
not considered very significant. Additional planes were tried to 
distinguish between other constants like Gy?, Cgg etc but without 
any success. It should be pointed out further that if the diffe- 

! t 

rence between two constants say C^g/^ ^88 etc. is not large com- 
pared to the inaccuracy (8^ for these constants) involved in the 
measurements, then they would not be detected and the results can 
be interpreted in terms of the classical theory. 

t f 

Evidently therefore, the difference between and Cgg and 
such differences in other constants ( constants involved 
should be determinable from individual K-values independently) 
can be employed to differentiate between the two theories, the 
limitation being that the difference should be larger than 10%, 
if the X-ray method is employed. Due to limitations of the 
accuracy of the method employed, it is not possible, on the basis 
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of our experimental results to make any assertion regarding the 
differences between the classical and new theory of elasticity, 
particularly for the case which has been investigated. 

5.1.3 Comparison of the experimental methods: 

5 . 1 . 3 . 1 Photographic method: 

In the present investigation only one relp (4 0 0) is studied 

by this technique .. Five misset positions on one side of the relp 

have been observed to give the elastic constant Gr-r-. The results 

55 

are in good agreement with those obtained by the diffractometer 

technique where equally small crystal for the same relp (4 0 0) 

i ? P 

was employed (e.g. C 55 - 3 .49(5) X 10 dynes/cm by photographic 
method, C^^=5.4S(3)x 10^^ dynes/cm^ by diffractometer). Each 
exposure for diffuse spot ran around 30 to 50 brs. The accuracy 
is limited due to the various other causes described earlier 
(Sec. 4 . 1 ). The photographic method takes very long time for 
the complete data collection compared to the diffractometer 
method. Results are given in Appendix I. 

5 . 1.3. 2 Diffractometer method: 

As pointed out in Sec. 4.2, small crystals of the size 
(ruio""^ cc) are generally used for the photographic methods and 
larger crystals of the size fu 5 mm x 5 mm '.rro (Rnnachandmn and 
Wooster 1951 ) customarily employed for the diffractometer tech- 
nique. In the present investigation a very small crystal of 
raiO“'^ cc has been successfully used not onlyki the nhotographic 
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method but also in the counter diffractometer method, The agree- 
ment between the results obtained by two methods with small crys- 
tals is considered good. As a resultjone of the limitations nor- 
mally associated with the counter techniq_ue is no longer consi- 
dered necessary. Further jcrystals varying in size by an order 
of magnitude led to essentially the same values of elastic constan 
when studied by diffractometer method (namely from ££*1 0 0 £[^qq 
and from I^[£l0 03 gQQ are same within the experimental limit). The 
additional advantage with such a small crjj’stal is that the diver- 
gence Correction is further reduced, 

5.1.4 Equi scattering contours around relps (4 0 0) and (6 00 ): 

Equiscatt ering contours are the surfaces surrounding each 
reciprocal lattice point over which diffuse intensity has a cons- 
tant value. These may also be called surfaces of equal diffu- 
sion. Jahn (1942) has shown how to calculate the approximate 
form of these surfaces. 

The method of drawing equiscattering contours for diffrac.- 
tometer and photographic technique is the same. These equi- 
scattering contours which are drawn for the given value of ■•the 
diffuse intensity becomes a particular case of K-surface calcu- 
lations (Sec. 5.1.5). 

The shape of equiscattering contours around the node (400) 
and (600) in the reciprocal, plane a*c* and a*b* respectively, 
is a typical form of Cassini’s oval family of curves, 'Consider 
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the equiscattering contour for node (40 0) for =0.77 x 10" 

in the reciprocal plane a c (Mg. 5.1a). The curve is symme- 
trical about a* and c*. The extension of the contour along a* 
shows that amplitude of vibration of the longitudinal wave is 
higher than transverse wave. Eq.ui scattering contours around 
node (6 00) in reciprocal plane a*b* drawn for I^/Iq =3.683 x 10" 
is also similar in shape to that around (4 0 0) node (Mg. 5. lb). 
It is symmetrical about a* an b*. let us now investigate the 
nature of these two iso-diffusion contours in relation to the 
mathematical expressions derived from the thermal theories of 
diffuse scattering of X-rays. Prom equation (3.1)^, We have 
”2 _ o2 kT JT- ™ 2 ;^2 ^ll^»^»^l!hkl 

where is corrected for all factors including polarization 


factor 
Or, 

node of a given crystal). 


-2 ■ K[3u,v,w(]. , . 

q = ( a constant ) x — (Por any particular 

■^ d ' "^0 


plane; 


hence, 


Now for an equi scattering contour around (4 0 0) in a*c* 


V = 0, and I^/Iq = Const. (1=1, 1=0, N=0) 


— (Const.) X ^00 Const. (A ^■j'] 


vn 
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whe re , 


Ai] 




((CllU +C^^w )(0^^u +G^^w‘")- (C^ 


... (5.1) 


= f('u^, w^) 


or q. = 0.f(u‘^5 w^) where G is a constant (5.2) 

f(-u , w ) a function of u and w given by Eq. (5.1). So the 
wave vector for a point on the contour is guided by the relation 
(Eq. 5.2), Since the function 'f' is dependent only on the squares 
of the variables u, w, (u and w being the direction cosines of 
the direction of propagation of the thermal wave corresponding 
to wave vector q). The value of q does not depend on the signs 
of u and w and hence equiscattering contour would be symmetrical 
about a* and c* axes and in general it would be an ellipse with 
its ax:es along a* and c*. let us now consider the values of q^ 
for ]j^ 0 13 and {Jl 0 0^ directions. Erom relations (Eq, 5.1) and 
(Eq. 5.2) and tables 3.1, 3.2, 3.3, we have, 


^001 ~ *^o^st. X f(0, w^)=Const. x 


(•.* w = 1) 


or 

and 

or 


^001 ^ ^55 


Constant . . . 


(5.3) 


^100 ^o^st. X f(u^, 0) 


qioo ^ ^11 = Constant 


Const .x^ 


’11 


(• , * u = 1) 


(5.4) 


Erom Eig. 5.1a, it may be easily seen that 

‘^100 ^ ^01 
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This inequality leads us to the conclusion (from Eqs.5.3 
and 5.4) that the elastic constants for this particular crystal 
would he such as to satisfy the inequality, > ^ 11 * ^ glance 

at tables in Appendix I for (4 0 0) relp vrould indicate that the 
results derived from the analysis of the equiscattering surface 
around the node (4 0 0) in a*c* plane, are in complete agreement 
with the elastic constants and vice versa. 


Now consider the nature and the shape of the equiscatter- 
ing contour around the node (6 0 0) in a*b* plane. The closed 
curve is symmetrical about a* and b* axes. 


Proceeding as before it would be seen that 
q^ = Const. = c.f (u^, v^) 


where - 


f{ut y2) = 




(O^^u^+CggV^) (OggU^+CggT^) 

... (5.5) 

2 P 

Thus q depends on f(u ,v ) whose value is independent of the sign 
of u and v. Consequently there would be symmetry of the 'isodiffu- 
Sion contours in a b plane. Further from Bq. (5. 5) jib is seen 
that , 

r2 _ V 1 (u = 0, V = 1) (5.6) 


1510 = Co^st.x^ 


i_ 

66 


Similarly, 


•^2 

llOO 


Gonst.x? 


'11 


(5.7) 
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From relations (Eqs. 5.6 and 5.7)s we see the inequa3-ity > 

Iqig ’ leads us to the conclusion tha,t elastic sonstant of 

this particular crystal would he such as to satisfy the inequality 
^66 ^ ^11* glance at tables given in Appendix I for relp (6 0 0), 
once again shows, that the results derived from an analysis of the 
equiscattering surface around the node (6 0 0) in a*b* plane are in 
complete agreement with the elastic constants and vice versa. 

The sensitiveness of the form of the surface to the elastic 
properties of the crystal suggests that the study of the diffuse 
scattering may become an important method of investigating those 
crystoAs which have marked anisotropy of elastic properties such 
as the present crystal, whose elastic properties depend in more 
marked way on the directions than do those of a cubic crystal. 

This variation of elastic properties along the surface is des- 
cribed in detail in Sec. 5.1.5. 

5.1.5 K-surfaces: 

The thermal diffuse scattering from a crystal varies markedly 
with direction. Along any^line (rekha) through a given reciprocal ^ 
lattice point (relp) the directional factor is denoted by | 

^C^’'^»^1]hkl ( I^^achandran and Wooster 1951), the value 
of which is proportional to the direction cosines (u,v,w) of 
the reMia and hkl,the indices of the relp. 

The expression for K may be written (Eq. 1.19), 
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where L = h/(h^+k^+l^)”''^^ etc. 

and ^ + CggV^ + etc. 


£y inserting particular values of elastic constants and giving 
(n, v,w) all possible values, the whole range of K- values may be 
covered. In the present investigation, the number of non-zero 
elastic constants is 9, as the symmetry is orthorhombic and 
expression for K contains correspondingly more terms. Although 
h, k, 1 may have any integral values, in practice very few 


values suffice for the study and representation of ela.stic pro- 
pert ies. JbJ? orthorhombic crystals O'] and '£00i] 

cover most of the ' requirements. 


5. 1.5.1 Construction of K- surfaces: 

For an orthorhombic crystal, we have studied K-surface 

for the relps (6 00 ) and (0 06 ) projected on the [^0 0 lj and 

[]0 1 03 respectively. These values of h,k,l and elastic constants 

are inserted into the formiuLa for K[^u,v,w[] evaluation 

is ca.rried out on IBM 7044 computer for large number of directions 

defined by the direction cosines (u,v,w). The values of K are 

plotted (3 times the actual scale) on a stereogram, small circles 

(constant 6-values) separated by 2° intervals and grea.t circles 

(constant if -values) also having a separation of 2° from one 

2 -1 

another are used. The dimensions of K are in cm -dyne , same 

as those of elastic moduli. K-contours are therefore numbered 
-15 

in units of 10 . 
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There are a number of applications of K- surfaces given 
in detail "by Prasad and Wooster (1956). The present investiga- 
tion demonstrates how K— surfaces are useful in showing which 
points on a given density contour are least likely to he affected 
h^r the second order correction (carried out in a similar way as 
the divergence correction). 

The density contour (solid curve) for first order thermal 
diffuse scattering for (6 00) relp projected on [^0 0 I'Q (5'ig. 
5.2a) shows a large K-value along Ql 0 C5~] and relatively very 
small K-value, along jjD 1 0~| • The K-value in any given . direction 
is proportional to the reciprocal of the corresponding elastic 
constants hence a large K-value (Pig. 5.2a) along {^l 0 0^ means 
a small value of elastic constant (C^^) and vice versa. The 
same is true for small K-value along [j) 10]} i.e. for (Cgg),this 
also shows that the thermal wave vector l-ioo^^OlO* Same conclu- 
sion was drawn from the equiscattering contour for the node 
(6 00) (Sec. 5.1.4 and Pig. 5.1h), 

The dotted line represents the K-surface (Pig, 5.2a) cor- 
rected for second order correction ( Ramachandran and Wooster 1951, 
Lucas 1968, 1969) for relp ( 6 00 ). Prom the contours (solid 
and dotted) over the whole stereogram, one can easily see the 
rekha K []1 OO^gQO affected most while K[}0 1 0]]5 qq ts 
c.ffected i-he least and the intermediate K— values (K|]]u, v,w]] |^^) 
are affected correspondingly. 



FlG.5-2(a)STERE0GRAPHlC PROJECTION OF THE K-SURFACE FOR 
POTASSIUM NI0BATEiRELP{6 0 0) PROJECTED ONfOOl] 
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Thus a large change is observed in a K-value along j^l 0 
after applying a second order correction (dotted curve), meaning 
thereby that the second order correction for is largest and 

least for Cgg. This is in excellent agreement with the corrected 
and uncorrected values and Cgg. Por convenience only a single 
contour is drawn with and without the second order correction. 
Similarly the contour projected on the ( 006 ) 

(Pig. 5 . 2 b) shows that is affected most by second order correc- 
tion fl/K COOOQQg = and the least (055=1/1^1 0 O^qq^) 

Thus the na,ture and shape of the isodiffasion contours around (4 0 0 
and (600) relps, and K-surfaces for relps (6 0 0) and (00 6) are in 
complete agreement with the theory of diffuse scattering of X-rays. 
The information given by the equiscattering and K-surface -studies 
regarding the relative values of the elastic constants is consis- 
tent with the values of the elastic constants of Potassium Niobate 
determined by X-ray diffuse scattering. 

5.2 Discussion 

5 . 2.1 Accuracy of the results: 

Accuracy of the results of the elastic constants, deter- 
mined by the X-ray diffuse scattering method depends upon the 
following factors, 

1. Choice of the reflecting plane (relp) and the rekha. 

2. Various corrections 2(i) absorption (ii) divergence 
(iii) skew, (iv) polarization (v) second and higher order, 

(vi) Mosaic, (vii) white radiation and (viii) general. 
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A * ’ 

FIG -5.2(5) STEREOGRAPHIC. PROJECTION OF THE K-SURFACE FOR 
POTASSIUM NiOBATEj HELP (0 0 6) PROJECTED ONC)I(0 


118 


scattering. 

5. Absolute intensity (^q) nieasurements. 

4. Cochran’s optical mode contribution. 

I'he first tv/o factors affect the accuracy of an indiV'idual 
constant in different v’sys, on the other hanr' the 3 ?lc-; ors, 3 and 4 
affect all the elastic constants to the nearly same extent and 
hence will be discussed latter. 

Constants C 22 aJ^d are determined from the rekhas 

^®spectively. These constants 

involve quite a large error pj 7% . The v;ave vectors for them 
are relatively large i.e. measurement is made away from recipro- 
cal lattice points, thereby making second order correction larger 
re 8% even if the higher orders are neglected. One could in princi- 
ple make observations close to the relp but the divergence correc- 
tions (particularly i divergence) increases rapidly nearer the 
relp and becomes uncertain. There are limitations due to little 
mosaic spread also, alternative way to reduce the second order 
correction is to obtain the elastic constants 0^^ etc, from relps 
of lower indices i.e. at smaller 20. This obviously increases 
skew correction to ru10% or so. In fact obtained from (400) 
relp does not show any appreciable difference from the one obtained 
from (6 0 0) relp. Higher the plane index higher is the diffuse 
intensity which improves the accuracy, and the structure factor 
value for these planes are also quite high. Therefore diffuse 
intensity can be measured with better accuracy. Divergence 
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correction has been reduced by using a small crystal rd 10“^ c. c. 
and small slits for beam and detector collimators. The absorption 
factor is the most serious source of error affecting the accuracy 
of the results in the X-ray work. This is reduced by using a well 
shaped crystal of suitable size, further the linear absorption 
coefficient of this crystal is very carefully determined (Appen- 
dix II I) and is accurate within 5 %. Errors due to mosaic spread, 
white radiations and optical vibrations etc. are corrected by 
taking mean of the intensity readings at points equidistant from 
the given reciprocal lattice point. All errors are further re- 
duced by studying the rekhas on either side of the Bragg reflec- 
tion and taking the average value of elastic constants. Therefore 
the accuracy of the C22 and is 7 %. The constants 

Cgg are obtained from the rekhas [)^l]Q^Q,[ino];QQg, [oio] res- 
pectively, Here the thermal wave vectors involved in the calcu- ; 
lations are very small compared to the case of C22 

Thus, second order correction is much reduced ( 2 %) while other 

corrections are nearly the same as for C22J Therefore 

the accuracy of the *^55 ®^d Ggg is 5 % • 

Constants 0^2* ^-13 ^23 obtained from the rekhas, | 

U i Hi: ^ ± ^lo06 f2^040 : 

respectively, for example, the constant determined as , 

follows , 
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similar is the case of 0 ^^ and The accuracy of these cons- 

‘tants is obviously low as they depend upon the values of the 
independent constants. Second order correction is q.uite appre- 
ciable in these cases ('^5%). ill other corrections remain 
practically the same as for etc. It may be noted that the 
accuracy of *^' 13 ? *^23 ^^.inly depends upon the accuracy with 

which the K-values for the relps are known, iny small error in 
the K- value introduces a large error in these elastic constants. 

An idea of the accuracy of K-values can be obtained by comparing 
the experimentally detennined K-values with those derived from 
the elastic constants, for example, the mean of the experimentally 
observed value K[7+ , O^^qq relp is given as 0.986 x 10“"^' 

cm^/dyne (Appendix I) (after the second order correction). 


kE 


+ 


f2’ 


1 

r2’ 


0] 


2(^22 ^ 66 ) 

1 ^^ 22'^^66 ^ ■^^ 22 ^ 66 “ 2'^1 2%6 



The value of the right hand side of this expression can be found 
by substituting the evaluated values of the elastic constants 
(Table 5.1). 


Thus, 

^ C± ^ » -t = 0*965 X 10“^^ cmVdyne 

The observed and calculated K-values agree well within the 
estimated error ( /^ 3 to 4 % )• The remaining two factors 
which reduce apcuracy of the elastic constants may now be 
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discussed. The accuracy of the determination of the elastic cons- 
tants from the study of diffuse reflections mainly depends on the 
accuracy v/ith which the ratio I^/l^ can he measured (in both dif- 
fractometer and photographic methods) separately in addition to 
the accuracy of the structure factor values. The measurement of 
the direct heai.o intensity is a very difficult task on the dif- 
fractometer owing to its very high value (rvlO^ times the diffuse 
intensity Therefore some intermediate standards are chosen, 

in the present investigation Compton intensity I , fix>m a dia- 
mond, which is of the order of diffuse intensity, is measured 
and from that the direct beam intensity is determined. The 
method is discussed in detail in the Appendix II. Such a stan- 
dardization of the incident intensity can be made with an accuracy 
of a few percent (ro3%). The accuracy of intensity measurements 
by photographic photometry is limited due to several factors 
(Sec. 4,^.). But in the present investigation I^ and 1^ have been 
measured simultaneously on the same photograph; The blackening 
for the diffuse spot and direct beam impression is adjusted or 
rather matched with the help of suitable numbers of Zr filters 
placed in the path of the direct beam. For measuring the inten- 
sity from the blackening of the photographic emolsion of the film, 
standard wedges were developed (simultaneously) with each photo- 
graph. The piece on which the standard wedge was recorded wan 
cut from the same sheet as was used for the corresponding diffuse 
intensity records (or monochromatic pictures). So it is believed 
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■fchat the accuracy in the measurement of is within 3% provided 

the liiiuar absorption coefficient of zirconium foil is absolutely 
correct. iQ.lowin;j for other errors like structure factor, optical 
vibration etc. , accuracy in the photographic method for the cons- 
tant is + 8%. 

Lastly we shall describe the deviation in the values of 
elastic constants due to the low frequency transverse optical 
phonon mode or socalled Cochran’s mode responsible for the ferro- 
electric property in the crystal. In the Sec. 1.4, a brief sum- 
mary of Cochran’s theory has been presented. The theory predicts 
the existence of a low frequency mode which goes to nearly zero 
as the critical temperature is approached in a given crystal. 

Such a low frequency optical mode would obviously contribute to 
the diffuse intensity besides the usual acoustic modes and hence 
the values of the elastic constants determined on the assumption 
of negligible optical contribution, would deviate from their true 
values. The recent study of Potassium Niobate by Perry et al 
(1969) by laser Raman spectroscopy does not show any such low 
frequency phonon mode in the room temperature orthorhombic phase. 
Other data such as the dielectric dispersion etc. which can throw 
light on the existence of the Cochran mode are not a,vailable. 
Further, the transition from ferroelectric to paraelectrio phase 
in the case of KlbO^ takes place near 435°C, while elastic cons- 
tants have been obtained at room temperature in the present study. 
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Thu ircQ.uuii cjf ira net abnorinally low al "this "t ©nip e nature (Penny 
et, al 1969). If one malces a veiy crude approximation that the 
cr;/atal ic cubic and follow the procedure outlined by Bell and 
Rupprecht (1963) for SrliO^, it is found that the deviation in 
tne elastic constants due to low frequency (Cochran's) mode is 
not more than 2 to 3%. Therefore taking Into account all the 
factors v/hich in fir ence the ■■ ccu:;tc'' of the elastic constants, it x 
is found that C^^, 022:1 C^^ can be measured with an accuracy of 
+ 10%, C^^ and Cgg with an accuracy of + while the, 

denendent cqnrtantr- ^•^ 2 ’ '"■^23 -'^13 ^ accurary of + \ 2 % , 

5.2.2 On the new theory of elasticity: 

The new theory has been discussed in Sec. 1.1.2. In Sec. 

5.1.2 attempts were made to determine the elastic constants on 
the basis of the new theory but it was observed that the experi- 
mental errors in the present method were larger than the diffe- 
rence between elastic constants based on the old and the new 
theories. In recent years, there have been few experiments per- 
formed to verify the new theory but it is surprising that even 
in the simple cubic crystals e.g. laBrO^, the difference in the 
elastic constants obtained by new and old theory is hardly 0.1% 
(Eadha and Rajgopal, 1968) which is well within the limit of 
experimental error of the method employed (pulse technique). 
Unequivocal results justifying the validity of the new theory 
have yet to appear. ■ 
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5*2«^ R&3.' '-ioi'iship oi" diffuse sccirl/Sriiig ■vri'i It oi’l)6X’ ’oi'op6!c4'v©.s 

liscc 'ohernal expansion, structure etc, 

A qualitative correlation between thermal expansion and 
diffuse sC'^ttering has been discussed in detail in Sec, 5*2,4» 

A brief mention of it may be made here. The diffuse scattering 
is entirely due to the amplitude of vibrations of the sfcoms, the 
larger the amplitude the greater is the observed diffuse intensity, 
Siiailarly thermal expansion is normally a function of lattice 
vibrations. Large theimal expansion coefficient along certain 
directions indicates a loose binding between the atoms i.e, the 
existence of wca3s; atomic forces or in other words, large amplitude 
of vibrations. This is what has been observed in study of diffuse 
scattering also, large thermal expansion is observed along 1 0'^ 
Similarly diffuse intensity is also having a comparatively large 
spread around (OkO) reciprocal point, lo studies on the thermal 
expansion of this single crystals are available to date to make 
any definite conclusion in the establishment of such relation and 
extending it to the individual values of the elastic constants. 
Among the elastic constants, it is worth mentioning that the 
direction along which the constants and are 

measured, undergoes a minimum change v/ith the temperature, as it 
can, be seen from the thermal expansion studies (Shirane et al 
1954) of the polycrystalline sample along this direction. Our 
studies of equiscattering contours and K-surfaces give ■ infor- 
mation which is consistent with the experimental values of the 
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elastic constants obtained in the present study. The Debye tem- 
perature of Potassium Niobate has been estimated, from the room 
temperature elastic constants, to be 302°K. The details of 


these computations are given in Appendix IT. The ratio of 


0.23. This value may be compared with 


that for iso structural BaTiO^ ~ = 0*21 (Tolger 


1952, Bhide et al 1965). 


5.2.4 The relation between the structure of ferroelectric 

Potassium Niobate: its thermal expansion and the X-ray 
diffuse scattering; 

In the present investigation an atten^jt has been made to 
investigate a relation between the structure of the Potassium 
Niobate and its thermal expansion behaviour in the orthorhombic 
phase, it is also shown that this relation is in full agreement 
with the observations of the thermal diffuse scattering. 

The graph (reproduced from Shirane et al 1954) ilg. 5.3 : 

shows the variation of lattice parameters with the rising tem- 
peratures, : 

On looking carefully into the graph (Pig. 5.3) of thermal | 

I 

expansion, a contraction is observed in the lattice parameters 
along a and c (a graph is plotted for Monoclinic parameters 
eQ_uivalent to given orthorhombic cell) sin 

2a’' cos p/2, b (symbols have their usual meaning) nhile 



A 



FIG* 5*3* Lattice paraJricter;; of as a fiinction 

of temperature. ‘'Ehir'-aie et al 1954) 



127 


there is considerably large expansion of the lattice parameter 
along the b direction, which is smallest (ro4°A) parameter, The 
rate of change of lattice parameters with temperature in the 
orthorhombic phase and the linear thermal, expansion coefficients, 
calculated from the data of Shirane et al (1954) are given here, 

= -0.944 X 10“^ i/°C, = -1.650 X 10"^ /°C 

|| = -0.222 X 10*"^ i/°C, = -3.896 x 10“® /°C 

II = 7.111 X 10"^ i/°C, = +1.791 X 10"^ /°C 

It follows from above that the linear thermal expansion 
coefficients are slightly negative for a and c, while it is mar- 
kedly positive for b-axis. The small thermal expansion coeffi- 
cients in (hOO) and (0 01) direction means that the binding 
forces along these directions are very strong, while the linear 
thermal expansion coefficient in the (OkO) direction which is 
large, shows that the atoms in this direction are loosely 
bound compared to other directions. We get this qualitative in- 
formation from the thermal expansion studies of a polycrystalline 
sample of KFbO^. The thermal expansion is highly anisotropic and 
its direct relation to the unit cell dimensions and bond directions 
seems quite difficult (Krishnan et al 1969). The directions of 
Tnlm*rrTOTn expansion (h 0 o) and (0 0 1) seem to have a close relation 
with dielectric anomaly (Wood 1951), while the direction of maxi- 
mum expansion suggests a large displacement of atoms fcom their 
equilibrium positions. 
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One can see "that a very careful study on the single crys- 
tal of this compound may throw light on the direction of Tna-yiTmiTn 
dielectric anamoly and mechanism of spontaneous polarization 
(Krishan et al 1969). 

Now the diffuse scattering studies on the single crystal 
of the same compound (KNhO^) supports the above observations, 
let us take the relps (4 0 0), (0 0 4), and (0 4 0). Their struc- 
ture factor is not too different (P^qq = 72.73; = 71.92; 

l'o 4 o “ 62.83). In the Pig. 5.4 we have observed variation of 
the diffuse flux or reflecting power for these relps vs. the 
misset from their respective Bragg positions on either side. It 
is clear from the Pig. 5.4 that the diffuse flux or power (04 0) 
relp is well spread over 5° to 6° on either side of the Bragg 
position for (0 4 0), while the diffuse flux for (4 0 0) and (0 0 4) 
relps is hardly observable after 3° on either side of their res- 
pective Bragg positions. 

As we know that the thermal diffuse scattering is entirely 
due to atomic vibrations, more the amplitude of vibration, more 
is the thermal diffuse scattering, therefore for (0 4 0) relp we 
see that atoms in this direction have relatively large amplitude 
of vibrations (displacements) than atoms along (4 0 0) and (00 4). 
Even the Bragg intensity for (0 4 0) relp is large compared to 

(4 0 0) and (0 0 4) relps ( 15 ( 040 ) > ^B(400) ^ ^B(004) 

shows that atoms in (0 4 0) having large amplitude of vibrations 




HO* 5,4* Variation of diffnae i|it#n«ity on altber eidt 
of tb« Srt«« foaltioni for C0 4 0)» (4 0- 0) 
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or birdinf: forces betv/een them are weal«ythar the atoms in (4 0 0), 
Similarly atoms in the (4 0 0) have slightly larger amplitude of 
vibrations than atoms in (0 0 4). The thermal expansion can be 
defined as follows. 

V/hen a solid is heated^ the atoms in the solid begin to 
execute vibrations about their equilibrium positions. The ampli- 
tude of these vibrations increase with the t eii 5 )e rature , If the 
atoms were executing purely harmonic vibrations the mean positions 
of the atom would not alter with temperature; hence solid would 
not expand. It is the anhaimonicity of these atomic vibrations 
v/hich causes a displacement of the mean positions of the atoms 
and thus brings about an expansion of the crystal. When the mean 
position of rest of an atom gets displaced, the forces acting on 
the atom tending to restore it to equilibrium position as it vib- 
rates, also alter and may cause a change in lattice frequencies. 

So expansion in the orthorhombic phase, means that since 
atoms along (0 4 0) direction have large amplitude of vibrations 
it should show large diffuse intensity, as is shown in the 
Pig. 5.4. Very recently some observations of this type are 
reported by Lambert, Gomes and Guinier (l968, 1969) and they 
found a disordered structure along (OkO) direction in the 
orthorhombic phase. 
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From tlio aniostropy observed in the tbermal expansion 
studies in tho orthorhombic phase of the polyciystalline sample 
of Kl'JbO^, it is very difficiilt to establish a quantitative rela- 
tionship with the observed elastic constants obtained from the 
diffuse scattering studies of a single crystal. Constants 
(longitudin;al component) and (transverse component) both 
having highest value along the respective directions and corres- 
ponding thcimal expansion in these directions is minimum. It 
shows that these constants are relatively temperature independent. 

Thus our studies of thermal expansion and the diffuse 
scattering ai-e in complete agreement, since they draw the same 
conclusion. 
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iPPENDIX I 


Calculations: Tables, Graphs and Plates. 

Table A- 1(a) 

CRYSTAIi DATA ON POTASSIIM NIOBATE 

Chemical Formula EUbO, 

3 

Crystal Class Orthorhombic |^mm2 


AXIAL PARAMETERS OF lOIbO^ 


a(i) 

b{i) 

c(i) 

Reference 

5 . 702 +.O 1 ' 

3. 984+. 01 

5.739+.01 

Wood ('1951) 

5.6946 

3.9714 

5. 7205 

Yousden (1951) 

T-* 

0 

> 

+1 

0 

• 

+1 

±•01% ■- 


5.6974 

3.9711 

5.7225 

Shirane et al (1954) 

5.697 

5.971 

5.721 . 

Katz and Megaw (1967)^ 

5.697 

3.970(9) 

5.721 

Present study 

Space { 

group 

Bmm2 

Number 

of molecules 

per unit cell = 2 


Yolume 

of the unit cell = 129.425 


Density 




Measured density 

= 4.62 gm/cc 


linear 

absoiption coefficient for 

MoKa 

Ae 

(calculated) 

= 58.207 cm"'*' 



(observed by diffractometer method) = 55.954 cm' 
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Description of the Tables A~2 to A-15 (Diffractometer 

Ohservations) 


Ccluran 

No. Description 

I. The angle of incident heam with certain given direction. 

II. Observed diffuse Intensity (units counts/sec. • or G.P.S.) 
Corrected for Absorption Jlosaic and white radiation 
corrections are applied, wherever necessary. 


III. 


Correction factor = 


Skew factor 
Polarization' factor 


X Divergence facte 


IV. Corrected intensity I^ (C.P.S.) due to 1st order. 

V. '.aucre of the wavr-lenoth ( '^) 'Of- "i b e el 0 . st 5 c ' wave 
relevant to the specific direction. 


(General Pormula: Symbols and Calculations, Diffractometer 
Observations), 


o 1 1 q 


■O 


hkl •* 


(3.1) 


or 


I jX 4 




The slope I^/^t obtained from I^ vs. gi^aph plotted 
for each rekha, after subtracting the background. * 
Therefore, 


2 

^ _ ^t » 
" K[:u,v,w_l^ “ * 


X (Const.) 


A-1 
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Where G 


liT fL^r,2 




pj X 


2 


G dependn upon structure factor (l^),rel Tector (i), solid 
angle (/I) and the volume (V) of the crystal. JU. ,the linear 
absorption coefficient end [>' the room temperature at vdiich 
observations are made, are both constant throughout the study. 


One gets the value of an elastic constant from 
knowing, G for a given rekha of a given plane. 

Below some relevant data required for planes studied ere 
given . 




Table 

A-l(h) 




Plane 

studied 

Axis of 
rotation 

Structure 
factor Pg, 

Rel. vec- 
tor I 
x10-f 
cm”*^ 

- Solid 
angle 
Jl , 
xIO^ 

Volume 
of the 
ciystal 

T X 105cc 

Constant 

11 

X 10~^ ' 

0 4 0 

1 00 

62.83 

1.0072 

9.797- 

1.838 

336.1 

6 0 0 

00 1 

41.44 

1.053 

9.797 

1.838 

159.8 

006 

0 1 0 

38.77 

1.049 

9.797 

1.750 

138.70 

4 0 0 

0 1 0 

72.73 

0.702 

1.382 

0.178 

3.323 

0 0 4 

0 1 0 

71.92 

0.699 

1.382 

0.178 

3.22 
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Table A-2(a5b) 

Rekha Loo"] 0^0 
(a) 


Anglo 0 " 
incider/l 
beam wifn 

rooii 

direcM^ion 

oVrr/'dir- 

fix?:e inten- 
sity correc- 
ted for 
Absorption 

X 10“^ 

Gorrec. 

tion 

Factor 

S.D/P 

Corrected*'' 
Intensity 
(due to ' 

1st order) 

X 10“^ 


19.89 

3.319 

1.178 

3.910 

2.889 

19.79 

2.547 

1.178 

3.001 

2.328 

19.69 

2.175 

1.185 

2.578 

1.980 

19.49 

1.647 

1. 169 

1.926 

1.482 

18.99 

0.997 

1.177 

1.173 

0.841 



Rekba 

C^^3o40 


(b) 





21.87, 

3.388 

1.225 

4.150 

3.969 

21.99 

2.561 

1.152 

2.950 

3.080 

22.09 

2.361 

1.098 

2.592 

2.575 

22. 19 

2.248 

1.030 

2.315 

2.209 

22.29 

2.246 

0.935 

2. 100 

1.827 

22.39 

2. 127 

0.889 

1.891 

1.580 

22.99 

1.208 

0.845 

1.021 

0.770 


n 

Mean value of the slope Ifl/A+ (less background) 

“51ar“ 


from the graphs 

10.0 

9.0 

Mean value of C^^(lst order) x 10”"*^ dynes/cm^ 

3.361 

3.734 

Mean value of 0.., corrected for second order 

X 10“ dynes/cm‘^ 

Average value of 0^^ = 3.5(58) x lo”^^ dynes/cm^ 
from Tables 2(a) and 2(b). 



3.372 

3.744 
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Table A-3(a, b) 


Rekha L0T0 ]q^q 


ia) 

Anglo 0 / 
incioG'!'.! 

team vvi't/i'i 

fooiJ 

directi -n 

Obc. dif- 
fuse 

intensity 
corrected 
for Absorp- 

tioii X lO”^' 




Correc- 

tion 

factor 

S.D/P 

Corrected 
Intensity 
(due to 

1 st order) 

X 10 “'^ 

(== 5 )^' 

X 10“^ 

21.87 

3.298 

1.226 

4.044 

0.600 

21.99 

2.212 

1.153 

2.550 

0.462 

22.09 

2. 183 

1.100 

2.401 

0.381 

22.19 

2.077 

1.055 

2.150 

0.315 

22. 29 

2.050 

0.951 

1.950 

0.272 

22.39 

1.9B8 

0.890 

1.770 

0.233 

(b) 


Rekha H^IOloiO 


19.97 

3.756 

1.230 

4.620 

0.462 

19.89 

3.178 

1.180 

3.751 

0.400 

19.79 

2.413 

1.181 

2.850 

0.333’ 

19.69 

2.016 

1.191 

2.401 

0.298 

19.49 

1.622 

1.170 . 

1.900 

0,210 

18.99 

0.865 

1.182 

1.023 

0. 132 

Mean value 

of the slope 

2 

1 a/\ 4 - (less background) 

"'“■“T(b)" 3(ir“ 

from graph 

G 

U. ' ' Kt 


73.0 66.0 

Mean value 

Of G 22 (1st order) x10“^^ 

0 

dyne s/ cm 

0.460 0.509 


Mean value of Q 22 corrected for second order 

X 10’*"^^dynes/cm2 0.492 0,540 

Average value of Ogg = 0.5(16) x lo”'^ d 3 ?nes/cm^ 
from Tables 5(a) and 3(b). 



01 
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Table A-4(a, b) 





Rekha [o, - 

" T2’ T^'^040 


Angle of 
incident 
beam with 

r 001 ] 

direction 

Obr. dif- 
fuse 

intensity 
corrected 
for Absorp- 
tion X 10~^ 

Corrected 
Correction Intensity 
Pactor I, (due 

to 1st 
order) 

X 10“'^ 

X 1 ( = 

X 10“^ 

21.87 

3.459 

1.227 

4.244 

1.984 

21.99 

2.610 

1.150 

3.002 

1.540 

22.09 

2.343 

1.089 

2.552 

1.289 

22. 19 

2.329 

1.031 

2.401 

1.104 

22.29 

2.326 

0.935 

2.175 

0.913 

22.39 

2, 100 

0.900 

1.890 

0.790 

22.99 

1. 176 

0.850 

1.000 

0.385 

(b) 


Rekha [p, 

1 1-1 

Y2» " Y2 J 040 


19.89 

3.241 

1 . 176 

3.812 

1.444 

19.79 

2.472 

1.179 

2.915 

1.164 

19.69 

2. 161 

1.180 

2.550 

0.991 

19.49 

1.625 

1.172 

1.905 

0.741 

18.99 

0.979 

1.180 

1.155 

0.421 





4(a) 4(b) 

Mean value 

of the slope la/Xt^^®®® 

background) 


from the graphs 


= 

22.000 18.50 


Mean value of I/kJO, - ^2^040 

l/K[p, y ^5 - Y^^040 order)x10“^^d3raes/cm^ = 1.527 1.816 

Mean value of l/Kp, and 

1/k[o, - =^3 corrected for second order 

X 10**^^ dynes/ cm2 = 1.669 1.759 

12 / 2 

Average value from Tabs. A— 4(a) and A-4(b) = 1,714x10 dynes/ can 

contd, . 
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Contd, . . 


From Tables (5.1), (3.2) and (3.3). 


^22^S3'^*^44^'^^33^44‘~^^23^44'‘^23 


K [O.+ ill. +Tl43nAn ^^°35 ®44^ 


i040 


or 


'23 


-C44 ± f cf^4-G^^C^^+(C^^+C,,)(G 


44-^33^44-^^35^44^ ^ 22 ' 


2 

m: 


72**^ 72^040 


) 


(Negative sign before tbe surd is never applicable) 
by substituting the values of C 22 a^<3. mean value of 

1/Kf0, ± ^ order). 

One gets, 


Mean value of G^^ <1st order) x 10-12 dynes/cm^ = --1.854 
Average value of 02 ^ corrected for second order 
calculated from Tabs. 4(a) and 4(b) = - 1 . 9 (l 4 )x 1 o”*^d 5 nes/cm^ 
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Table 

A-5 (a,b) 


La). 


Rekha 

010] 600 


infio c.r 
incideji'b 
beam wi+.l" 
[010] 

direct 0>ii 

Ob.'^. dif- 
,fuco 

intenrity 
corrected 
for Absorp- 
tiOL , 

X 10~^ 

Correc- 

tion 

Factor 

S,D/P 

Corrected 

intensity 

(due 

to 1st 
order 

X 10”'^ 

/\t 

q. 

X 10“^ 

21.10 

3. 105 

1.063 

3.301 

3.906 

21.03 

2. .'92 

1.118 

2.786 

3.249 

20.96 

2. 147 

1.125 

2.416 

2.809 

20.90 

1.915 

1.163 

2.227 

2.525 

20.80 

1.573 

1.177 

1.852 

2. 116 

20.90 

0.834 

1.200 

1.001 

1.332 

(b) 


Rekha 



22.95 

2.383 

1.100 

2.621 

3.164 

23.00 

2.071 

1.038 

2. 150 

2.862 

23. 10 

1.759 

0.997 

1.754 

2.352 

23.20 

1 . 428 

1.012 

1.445 

1.936 

23.50 

0.909 

1.110 

1.000 

1.296 


Mean value of the slope (less backgroimd) 

from graphs 

5(a) 

- 8.0 

5(b) 

7.62 

*1 P P 

Mean value of Cgg (1st order) x10“ dynes/cm 

= 1.997 

2.096 

Mean value of Cgg corrected for second order 

X lO"""^^ dynes/ cm2 

= 2.02 

2,118 

Average value of Cgg = 2, 0(69) x lo''^ dynes/cm^ 
from Tables 5(a) and 




5(t). 
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Table A-6(a,b) 


(a) 


Rekba [_100^ 

600 


Angle 0: 
incident 
beam with 

roio] 

di 1*60 tic n 

Obfi. dif- 
fuse 

intensity 
corrected 
for Absorp- 
tion X lO"*"^ 

Correc- 

tion 

factor 

S.D/P 

Corrected 

Intensity 

(due)^ 

to 1st 
order) 

X 10“'^ 

(= 

X 10-5 

21.10 

2.328 

1.065 

3.012 

0.600 

21 . 05 

2. 139 

1.192 

2.550 

0.529 

20. 9o 

1.908 

1.127 

2. 150 

0.46 2 

20.90 

1.693 

1.165 

1.973 

0,400 

20.80 

1.276 

1.185 

1.512 

0.342 

20.70 

0.974 

1.213 

1.182 

0.289 



Rekha pool 

600 


(b) 





22.05 

3.092 

1.300 

4.020 

0.638 

22.90 

2.671 

1.125 

3.005 

0.576 

22.95 

2.380 

1.099 

2.616 

0.506 

23.00 

2.167 

1.039 

2.252 

0.473 

23.10 

1.738 

0.998 

1.735 

0.390 

23.20 

1.478 

1.015 

1.500 

0.32/h 


2 

Mean value of the slope I./'^+(less background) 
from graphs ^ 

Mean value of C^^(1st order)x10“''^dynes/cm^ 

6(a) 

49.5 

0.323 

6(b) 

47.5 

0.336 

Mean value of C^^ corrected for second order 

X 10""^^ dynes / cm^ 

0.347 

0.361 

Average value of « 0.3(54)3d0^^ dynes/ cm^ 

from Tables 6(a) and 6(b). 
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Table A-7( 

a,b) 



Xal 


Rekha Cps “ 

' W °^600 



Angle of 
incidcnf 
bean witl; 

r 010 ] 

direction 

Obc. dif- 
fiise 

intensity 
corrected 
j'or Absorp- 
tion X 10-4- 

Correc- 

tion 

factor 

S.D/P 

Corrected 

Intensity 

(due 

to 1st 
order) 

X, 10"^ 


(= 

X 10“^ 

22.90 

2.573 

1.128 

2.902 


1.755 

22.95 

2 . 36 2 

1.101 

2.601 


1.582 

25.00 

2.212 

1.040 

2.301 


1.437 

23.10 

1.747 

1.030 

1.800 


1.176 

23.20 

1 . 485 

1.010 

1.500 


0.968 

23.50 

0.913 

1. 150 

1.050 


0.648 

(b) 


Rekha [!• 

W °^600 



21.10 

2.934 

1.065 

3.125 


1.953 

21.03 

2.377 

1.115 

2.651 


1.625 

20.96 

2.066 

■ 1.125 

2.324 


1.404 

20.90 

1.853 

1. 160 

2. 150 


1.263 

20.80 

1.441 

1.180 

1.700 


1.058 ■ 

20.70 

1.101 

1.226 

1.350 


0.905 

20.50 

0.843 

1.210 

1.020 


0.666 






“TX^ 7 CR 


2 

Mean value of the slope I^/A+ (less 'background) 

from graphs ^ = I6.4 15*6 

lean value of ^"^600 

I/kQ- Y^* ^600 ^ ^ order)x10 ^ dynes/cm = 0*974 1*024 

Mean value of l/Kj^j^,- p’^BOO 

1/k;[- Tr^>nfi"»*3^nn corrected for second order 

X 10“'*^ a^es / cm^ 0.983 1.045 

Average value from tabs. 7(a) and 7(b)=1.0(14) x 10 dynes/cm^ 

contd* . 



U8 

Contd. . . 

Prom Tables (3.1), (3.2) and (3.3) 

V.li I 

or 

^12 “^ 66 i ^^ 66 '^^ 22 ^ 66 '^^^ 22 ^^ 66^^°11 “ ) 

K[± Y5» +72>?l600 

(Negative sign before the surd is never applicable) 

Substituting the values of ^■]-|»^22’%6 mean value 

of 1/K[+ Y^, 03 gQQ (due to 1st order). 

One gets, 

The mean value of order) x10“ dynes/cm = -1.069 

Average value of corrected for second order 
calculated from Tabs. 7(a) and 7(b) = -i.o( 58 )x 1 o'^^d 3 mes/cm^ 
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Table A-8(a,b) 


Rekha [100] 5 

(a) 


iBgl-;/ O': 
in c Hie ''ll 
be.'.'un wit.il 

Dool 

direotii,)':. 

Gbs. dif- 
fuse 

intensi'i'.y 
corrected 
for Absorp- 
tion X 10”^ 

Correc- 

tion 

Factor 

S,D/P 

Corrected 
Intensity /\ 

(due 

to 1st 
order) 

1 ( = 

q. 

X 10 




X lO”^ 


22.60 

1.564 

1.279 

2.000 

4.692 

22.73 

1 . 387 

1.255 

1.741 

4.225 

22.83 

1 . 092 

1.193 

1.303 

3.364 

22.93 

0.840 

1.244 

1.045 

2.782 

• 

CM 

0.669 

1.204 

0.805 

2.093 



Rekha 000] 006 


(b) 





21,00 

1.449 

1.314 

1.905 

4.624 

20.98 

1.053 

1.329 

1.400 

3.540 

20.93 

0.954 

1 . 366 

1.303 

3.192 

20.03 

0.816 

1.244 

1.015 

2.652 

20.43 

0.579 

1.079 

0.625 

1.369 





8(a) 8(b) 

Mean value of the slope 
from graphs 

IdA-t (less background) 

3.89 3.80 

lean value 

of 0^^(lst order)z10'“ 

dynes/ cm 

3.566 3,651 

Mean value of Gc-c corrected for second order 

- 1 9 9 

X 10 dynes / cm 

3.610 3.691 

Average value of Cj.^ = 

3.6(51) X 

10**^ d3mes/cm^ 


from Tables 8(a) and 8(b), 







152 


Table A-9(a,b) 


(a) 

Angl-; Oj" 
inoionni, 
beam v.'it'i: 

r 

direct 

Obs. d:^-. 

.t\j ;;.e 

i nte.'oity 
corrected 
for Absorp- 
ticc X 10 ~‘T 

Rekha [ 001 

Corrected 
Correc- Intensity 

tion I (due 

Factor ,, , 

S.D/P 

order) 

X 10 “^ 

XI (= 

q 

X 10 “^ 

22.69 

1 . 4 .; 6 

1.279 

1.850 

0.756 

22.73 

1.r?G9 

1.265 

1.605 

0.576 

22. B3 

1.040 

1.202 

1.250 

0.506 

22.93 

0.693 

1,226 

0.850 

0.421 

23. 13 

0.537 

1.210 

0.650 

0.225 



Rekha £ 001 ]oq 5 


(t) 





21,09 

1.381 

1.304 

1.801 

0.756 

20. 9 H 

1 . 032 

1.309 

1.351 

0.689 

20.93 

0.901 

1.360 

1.225 

0.552 

20.83 

0.598 

1.250 

0.748 

0.441 

20.43 

0.536 

1.085 

0.582 

0.324 


p 

Mean v^aue of the rjlope I, /X^ (less background) 
from gra.phs 

Mean ¥alue of (ist order)x10”^^dyes/cm^ 

Mean value of corrected for second order 
X 10”^^ dynes/cm2 

Average value of from. Tabs. 9(a) and 9(b) 


9(a) 9(b) 

- 23.5 22. 15 

= 0.590 0.626 

= 0.614 0.650 

=0.6(39)x10^^^. 
dyne s/cm^ 
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Table A-10(a,b) 


(a) 

Rekha 



Angli2 c.f 
iuciie-it 
beam v,l-i:h 

rwo,-] 

direc-”,\.M. 

ObG. dif- 
fuse 

intoaoity 
cori'-vTcted 
for Absorp- 
tion X 10"'^ 

Correc- 

tion 

factor 

S.D/P 

Corrected 

Intensity 

(due 

to 1st 
order) 

X lO"'^ 

1 (= ■ 
q 

X 10“ 

22.69 

1 . 526 

1.279 

1.952 

2.345 

22.73 

1.381 

1.260 

1.740 

2.112 

C\l 

CM 

1.088 

1.195 

1.300 

1.682 

22.93 

0,762 

1.254 

0.955 

1.391 

23. 15 

0.618 

1.214 

0.750 

1.046 

(b) 

Rekha 

1 

''^ 006 


21.09 

1.434 

1.310 

1.879 

2.311 

20, 9B 

1.034 

1.329 

1.374 

1.770 

20,93 

0.916 

1.367 

1.252 

1 .596 

20.8:) 

0.683 

1.245 

0.850 

1.326 

20.43 

0.540 

1.082 

0.584 

0.684 



10(a) 

10(b) 

o 

Mean value of the slope I^/>v^(less background) 



from graphs 

= 8.0 

8.0 

Mean value of 1/Kj}- snd 

1/K|^,0,~ ('•s'fc order)x10"^^dynes/cm^ 

Mean value of 1/K[1 ^ 

= 1.734 

1.734 

1/E;[ 0,- y?^ 006 °03n?ected for second order 

X lO"”^^ dynes/cm^ 

= 1.786 

1.786 


12 2 

ATCrage value from Tabs. 10(a) and 10(b) =1;7(B5)x:10 dynes/cm 



PrcE tables (3.1), (3.2) and (3.3), 



11 




)+C 


2(0 


11 


55 1 3 JJ_ 

+ C55-) 


2 

13 


^13 ■'^55 ~ ^55'^S3%5^^^33'^%5^*'°11 " ~”“™“ ]— 

^D: ir5»o»+ 7^1 006 

(Negative sign before surd is never applicable) 


Substituting the values of and the mean 

value of 1 /kI; + y 5 , 0 , +Y 1 I 0 O 6 (due to 1st order). 

One gets, 

1 P p 

The mean value of C^^(ist order)x10'' dynes/cm = -2.20 

Average value of G^^ corrected for 2nd order calculated 
from Tabs. 10(a) and 10(b) = -2.0(27)x10^^ dynes/ cm^. 
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Contd, . , 


Prom tables (3.1), (3.2) and (3.5), 


1/KCi ^,0,4- 


or 


0^3 = -C55 ± - 


K[+ ^||,o,+ Y^looe 

(Negative sign before surd is never applicable) 


Substituting tbe values of C^^,G^^,G^j and the mean 
value of 1/K[^+’^,0, +Yklio06 order). 

One gets, 

i P p 

The mean value of G^^(lst order)x10'' dynes/cm = -2,20 

Average value of corrected for 2nd order calculated 
from Tabs. 10(a) and 10(b) = -2.0(27)x10‘'^ dynes/cm^. 
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Table A-1l(a, 

b) 


(a) 


Rekha [TooJ^qq 


Angle o;*" 
incilen-i- 
■beaJ!- with 

r oc'* 1 

airection 

Obr. dif- 
fuse 

intensity- 
corrected 
for Absorp- 
-3 

tion X 10 

Correction 

factor 

S.D/P 

Corrected 

Intensity 

(due 

to 1st 
order) 

X 10“'^ 

* 

X 10”^ 

13.54 

C.635 

0.895 

0.568 

0.564 

13.44 

0.389 

0.990 

0.585 

0.441 

13.34 

0.303 

0.992 

0.501 

0.361 

13.24 

0.224 

1.006 

0.225 

0.298 

(b) 


Rekha [j OOj 

400 






15.34 

0.548 

1.103 

0.604 

0.529 

15.44 

0.392 

1.149 

0.450 

0.441 

15.54 

0.297 

1.185 

0.352 

0.324 

15.64 

0,225 

1.220 

0.275 

0.306 


Mean value of the slope I^/>^(less background) 
from graphs 

11(a) 

(b) 

1.0 

1.06 

lean value of C^^(lst order)x10“ dynes/ cm 

0.332 

0,314 

Oecond order 'correction is very small hence 
neglected. 

0.332 

0.314 

Average value of 0^^ =0.3(23) x 10^^ dynes/cm^ 
from Tab. 11(a) and 11(b). 





X 10 Id X 10 
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Table A-12(a,b) 
Rekha tpOll 4 oo 


( S') 


Angle of 
incident 
beam v/ith 

Tool ] 

direction 

Obs. dif- 
fuse 

intensity 
corrected 
for Absorp- 
tion X 10~5 

Correc- 

tion 

factor 

s.d/p 

Corrected 

Intensity 

(due 

to 1st 
order) 

X 10“^ 

X 10~^ 

13.54 

0.917 

0.892 

0.818 

8.190 

13.44 

0.612 

0.989 

0.605 

6.520 

13.34 

0.536 

0.990 

0.531 

5.220 

13.24 

0.445 

1.005 

0.447 

4.761 

12.44 

0.154 

1.012 

0. 156 

1.661 

(b) 


Rekha |p01 J" ^qq 


15.34 

0.755 

1.101 

0.831 

8.10 

15.44 

0.565 

1.148 

0.649 

6.52 

15.54 

0.452 

1.184 

0.535 

5.18 

15.64 

0.391 

1.219 

0.477 

4.49 

16.22 

0.149 

1.304 

0. 195 

2.07 


12(a) 12(1)) 

Mean value of the slope I,/-X^(less hackground) 


from graphs 0.094 0.098 

Mean value of (1st order)x10”'*^ dynes/cm^ 5.535 3.391 

Second order correction is negligibly small 

Average value of =3.4(6?) 'X 10 dynes/cm'^ 

from Tabs. 12(a) and 12(b). 

This value of 0^^ may be compared with those obtained from rekhas 
p 0 0^006 L'^^Olooe A':8(a, b)), where a comporptively 

I'ar'Cie ci“^/stal is used (¥ol. 1.750 x 10 ' cc,). 



01 X 




Table A-13(a,b) 
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Rekha £- 


(a) 

Angle of 
incident 
beam with 

rooi j 

direction 

Obs. dif- 
fuse 

intensity 
corrected 
for Absorp- 
tion X 10~ > 

Correc- 

tion 

Ractor 

S.D/P 

► 

Corrected 

Intensity 

(due 

to 1st 
order), 

X 10-5 

Si 

X 10“5 

-13.54 

0.860 

0.893 

0.768 

4.074 

13.44 

0.594 

0.987 

0.586 

3.261 

13.54 

0.505 

0.991 

0.501 

2.610 

13.24 

0.425 

1.004 

0.427 

2.381 

12.44 

0.139 

1.013 

0.141 

0.830 



Rekha » 

" TI3I4OO 


15.34 

0.728 

1. 102 

0.802 

4.051 

15.44 

0.531 

1.147 

0.609 

3.262 

15.54 

0.426 

1.184 

0. 505 

2.593 

15.64 

0.358 

1. 219 

0.437 

2.244 

16.22 

0.127 

1.307 

0. 166 

1.047 





13(a) 13(b) 

Mean value 

of the slope 

(less iDackgroiind) 


from Graphs 



= 0.18 0.19 

Mean value 

of 1/Z[1 

°’T2 3400 



7^^400 

order) 

dyaes/cm^ 

= 1.846 1.748 


Second order correction is negligibly small. 

Average value dErom Tabs, 13(a) and 13(b)= 1 . T(97)3:1o"*^dynes/cm^ 


This value may be compared with the value obtained your rekhas 

± f2 » + 7^Iloo6 ("tables A-10 (a,b) since they give the same 

constant 







X 10 


1-5H 

3 


O.R?- 



2 4 ~ '6 

,2 ,_~3 02 

A|XlC A 


p-- 


SECTION 1 TO [OlO] 



,2 .^-3 02 

At X 10 A 



165 


Table A- 15 
Rekba 


iingle of 
incident 
beam with 

r 100 3 

Erection 

Obs. dif- 
fuse 

intensity 
corrected 
for Absorp- 
tion X 1Cr> 

Correc- 

tion 

factor 

S.D/P 

Corrected 

Intensity 

(due 

to 1 st 
order)„ 

X 10 “^ 


(= Jpa' 

a 

-3 

X 10 

13.48 

0.848 

0.892 

0.756 


4.096 

13.38 

0.666 

0.922 

0.614 


3.342 

13.23 

0.553 

0.936 

0.518 


2.775 

13. 18 

0.439 

0.974 

0.428 


2.313 

12.38 

0.150 

0.982 

0.147 


1.601 

Mean value 

of the slope I^/ A ^ (less 

5 background) 

= 

0.18 


from the graph 

Mean value of I/kI^y^jO, - Y^jl] 004^ order) 

X 10 "^^ dynes/ cm^ = 1.7 ( 91 ) 

Second order correction is small hence neglected. 

'L'liif'^ value may be compared with the value of 1/k[I+ Oj + ^^^006 
which is 1.7(86) x 10^^ dynes/cm^, corrected for second order 
(table A- 10 (a,b)). 
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Appendix I Gontd. . 


Photographic Method 

Pew curves (Pig. 1-14 (a, h)) showing the variation of 
the intensity of diffuse reflection expressed as a fraction of 
incident intensity and uncorrected for polarization and skew 
corrections i.e. with the angle of diffraction 0 observed 

in the Lauo photographs taken with monochromatic MoKa (Plate III) 


radiation. 
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Tabic A -16 

(Photo graphic Fethod) 


[^0 1 oQ axis vertical and coincident with the ajcis of 
•the cejncra. ( Radius = 5.70 cm). 


Index of the plane under study 
Structure factor of the plane 
Thermal wave vector relevant to the table 
Temperature at which study is made 
(Room temperature) 

’lavelcngth of the monochromatic MoKa 


(400) 

72.73 

’^1.400 

300°IC 


radiation 0.7107 i 

Volum.c of the crystal 0.996 x 10"'^' c.c. 


Rekha £0 0 1 ] 


Angle of 
incident 
beam with 
[0 0 1 ] 
direction 

2x10-3 

Angle of 
diffrac- 
tion 0 
in degrees 

^ X 10^ 

0 

Corr, 

factor 

1 X 10« 
■^0 

C^^xlO"^^ 
d3?nes/cm^ = 

Vk[ooi3^oc 

13.425 

6.006 

28.89 

2.579 

1.051 

2.711 

3.472 

13.337 

5.476 

28.88 

2.315 

1.045 

2.418 

3.549 

13.275 

4.900 

28.86 

2.107 

1.041 

2. 194 

3.500 ! 

13.235 

4.692 

28. 88 

2.039 

1.040 

2. 120 

3.469 


The intensity 7 ^/Iq is corrected for general scattering. 

12 

Average value of 0^^ from the table A-I 6 is 5.497 x 10 

dynes/cm^ 

12 

This value is in fair agreement with the value of 0^^ = 3.463 x 10 
dynes/cm^ obtained by diffractom.eter method from the rekha 
Qo 0 i"34oo (a>’b)) , on equally small crystal 

( Vol. 1.782 X 10“^ cc.). 
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Lt> j 



PlilE I 

n-F ■pntaqsium Uio’ba.’te about "b-axis 
Rotation photographs of Potassium « 

+ ,» Cpi-vstal reflsotea) HoKa radiation, 

A,^ MonoctiT'oiiS't ic (. y 7Pp t i 0 ’”^oo, 

volume of the crystal 1.782 x iu o . 

volume reflected) ,Cu£a radiation, 

h) Monochromatic (crys ^ 

volume of the crystal 0.99b x lu 


PLATE II 


Orientation pictures (i.e. Laue photographs using poly- 
chromatic X-ray beam from Mo anticathode) with j^O 1 0^ 
axT.s vertical and incident beam makingj 

a) 13,275°; (4 0 0) diffuse reflection on the eq.uatorial line, 

b) 13,337°; (4 0 0) diffuse reflection on the equatorial line^ 
with c* axis. 


Ca) 
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PLiTE III 


Laue photographs using monochromatic MoKa radiation 
with [jb 1 CQ axis vertical and incident beam making 

(Shown by ^■nowhta.<() 

a) 13.275^; (4 0 0) diffuse refLection^on the equatorial line^ 

b) 13^331^; (40 0) diffuse reflection^on the equatorial line^ 


with c* axis 
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iLPPMDIX II 


MEASUREMEOT? OE EIREGT BEAM lETElfSlIY (Iq) 

To determine fhe elastic constants, it is necesssxy to 
measure the absolute value of the incident flux. Ihe direct 
X-ray beam is too strong to be measured directly with a pro- 
portional or scintillation counter, as it may not possess a 
linear response at these high counting rates; it must be atte- 
nuated first so that the dead time of the counter is appreciably 
less than the mean time interval between the arrival of the X-ray 
q.uanta. To evaluate the intensity l^, the attenuation factor 
must be accurately known. Multiple foils can be used for atte- 
nuation, but difficulties are caused by the progressive harden- 
ing of the beam, due to preferential absorption of the softer 
components in the passage through successive foils. The multi- 
ple-foil technique is, therefore, most suitable for measurements 
with mono chro mat i zed radiation or with balanced filters. This 
technique is capable of establishing the intensity scale factor 
to within 3 % of its true value (Burbank, 1964-65) . Since the 
reduction factor is of the order 10^, small error in measurement 
of thiclQiesses of the foils or /I, the linear absorption coeffi- 
cient, may lead to large errors in the determination of 1^, the 
direct beam intensity (Wooster, 1962). 

Another method of determining 1^ (Buyers, 1964) consists 
of measuring the intensity scattered by a small block of paraffin 
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placed in the main beam, and using theoretical values for the 
absolute magnitude of this scattering (international tables for 
X-ray Crystallography). In an amorphous solid such as paraffin 
the scattering from each atom bears a random phase relation to 
that from other atoms; thus the total scattering is contributed 
by independent unmodified scattering and by modified (Compton) 
scattering, and both of these can be calculated from the immodi- 
fied and modified scattering intensities of the separate atoms. 
The method was abandoned since the purity of the chemical com- 
povsition of the wax was not known correctly. 

Another intermediate standard can be provided by the 
Intei'.sity of the Compton scattering from diamond (Prince and 
Woostei' 1953) or other materials containing elements of low 
atomic number (Paraffin, etc._ described above). The whole of 
the incident X-ray beam is allowed to fall on the diamond crys- 
tal set at any convenient angle (9^^ -I or 6220^ incident 

beam. The radiation of modified wavelength scattered in direc- 
tions well away from Bragg reflections (about 8 to 10° on either 
side) has an intensity which may be calculated from the funda- 
mental constants and may be compared directly with the thermal 
diffuse scattering under investigation. ¥e may define a Compton 
diffuse scattering power, I)Q,in an analogus manner to the first 
order thermal diffuse scattering power, i.e. as the ratio of 
the intensity scattered by the Compton effect per unit cell of 
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the crystal per irait solid angle to that scattered by a single 
Thomson electron mader the same conditions (Wooster, 1962). 
For diamond this Compton diffuse intensity, I)„, may be written 
(James, 1948), 

= 8(6-5: (11-1) 

where Z is a function, which descrihes the effect of the 
individual carbon atoms scattering incoherently. Compton and 
Allison ( 1935) give a table listing the values of E ff as a 
function of the angle of deviation, 29, of the X-rays. It is 
also given in International tables for X-ray Crystallography 
Vol. III. B is Breit-Dirac correction factor, the value of 
whiob is, 

TO . . 2hX Sin^e 

B _ 1 +-5^ -y- , 

where the symbols have their usual meaning. 


From this calculated Compton diffuse intensity, D^, we 
may determine the absolute intensity of the incident beam I^, 
using the measured intensity of the diffuse beam, I^jand a for- 
mula analogous to that of the diffuse flux, namely. 


ic “o-m 


(II-2) 


here is linear absorption coefficient of diamond for MoXa 
wavelength, and V is the volume of the unit cell, other symbols 
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are explained earlier. The experiment is performed with balan- 
ced (Zx-t) filters as described earlier (Sec, 4»4- 1.1 .) 


Description of Table AII-1 


Column Description 

1. The angle which incident X-ray beam mahes with 
diamond crystal. 

2. The value of 47t Sin e/XioKa* 

5. The value of 1/B^ (Breit- Dirac factor). 

4. Polarization factor. 

5. Observed Compton diffuse intensity, I . 

c 

6. Compton diffuse scattering flux, D . 

7. Direct beam intensity, 

Thus for a brass collimator pair the average value of 

O 

the direct beam intensity is 4.886 x 10 c.p.s. Similar experi- 
ivan done for G. E. Collimator pair, the average value being 


merit v 


5.26 X 10^ c.p.s. 


I‘or diamond, 


A, 


Y 


1.9375 cm"' 
4.50 X 10“^^cc 


-TL, = 9.797 X 10' 

brass 


XL 


G.E. 

e2 


1.382 X 10' 

7.935 X 10 


-3 

-3 

-26 / 1+cps^ 29 


( 
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Table AII-1 


6° 

4-TCSine 

XploKa 

1 

1+cos^20 

2 


^c 

1^x10"® 

24.375 

7.300 

0.966 

0.717 

67, 135 

38.640 

4.676 

25.375 

7.580 

0.963 

0,700 

65.615 

38.520 

4.870 

26.375 

7.858 

0.960 

0.683 

68.625 

38.784 

4.995 

27.375 

8.134 

0.958 

0.667 

67.610 

39.086 

4.999 


Average value of the direct beam intensity, I 

from this table (for Brass Collimator pair) = 4.886 x 10® c.p.s. 
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iPPMDIX III 


DBTBB'IITATION OB IIBBAR ABSORPTIOF COEBBICIBxa' OB 

mo^ 

(i) Photographic Method ( Chakra'b'brty 1958 rnd SrivaataTa i960); 

A small crystal was moimted on a goniometer head. After 
adjusting one of its set of parallel and opposite faces perpen- 
dicula.r to the incident 'beam (Monochromatic and without harmo- 
nics), the direct heam impression was recorded on a piece of 
film '.Big-AIII (b) . This impression has at its centre the sha- 
dow of the crystal and unabsorbed direct beam on either side 
of this shadow. This piece of film was developed simultane- 
ously with a standard wedge and both were microphot ometered 
under the same conditions,? Big. ATII( c) . The scanning spot of 
light wos kept quite small such that it was much smaller than 
the shadow of the crystal and also the unabsorbed direct beam 
on oithor side of the shadow. The microphotometer record ha,s 
a dip in the centre with shoulders on either side Big. .,^111 (a) . 
Thf; intensity corresponding to the shoulders will give and 
that corresponding to the dip will give Thickness of the 

crystal t traversed by undeviated beam can be accurately deter- 
mined by measuring the distance between the pair of the parallel 
facor: which were made perpendicular to the incident beam. 



*1 






(b) 



no. A,III (a) Direct 1}eam iirpression with the crystal 
tn the path. 

(h) Direct beam ismressioh without ' the cryitH. 
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Hence /^can He evaluated from the following relation, 

>4.= log (III-1) 

^ a 

The linear absorption coefficient of Potaesium Hiobate 
was thus determined as 50.2 cmr''.jXcan also be calculated from 
the relation, 



(I II- 2) 


where is mass absorption coefficient of individual 

atoms consituting the crystal (International tables) p^ their 
proportion and ^ the density of the crystal. The value thus 
calculated for the present case came to be 58.207 cm””^. Inci- 
dent! y it may be mentioned that the agreement between the 
observed and calculated value of linear absorption coefficient 
(of KhbO^) in the present investigation is fairly good compared 
to the data reported by other workers (Brill et al 1939 axid 
Srivastava i960 etc.) on other crystals. 


(ii) Diffractometer method; 

A thin crystal plate of KFbO^ with opposite faces parallel 
to each other is mounted on the goniometer head which is in turn 
mounted on the single crystal orienter. A reflection from the 
plate at very low angle, close to O°(0) is balanced with Zr-T 
pair O'- niters. The crystal plate with its opposite faces 
parallel in brought in the path of the direct beam such that 
these parallel faces are perpendicular to the direction of the 
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incident team from Mo target, operated at 35 EYP and 25 ^lA. The 
detector is also in line with beam collimator i.e. the position 
20 = 0° on the 20 circle. The intensity of the direct beam (I^) 
is reduced to after traversing the thickness of the plate, 
t cm. This is recorded with the balanced pair of filters (Zr~T)» 
Several observations of the time for 10,000 counts are taken for 
both the filters. 

The crystal plate is then removed from the path of the 
direct beam (Iq). Keeping the detector in the fixed position 
(20 = 0°), the direct beam intensity is measured with the same 
balanced pair. 

li?om the formula ($1. A III-1) 

I = I e“ 
a 0 

„ M-t ^ _ j 

or e 

/Ut = log 

a 

1 ^0 

In this way Jit , the linear absorption coefficient for 
KKTbO^ crystal, is found. Observations are carried out on a 
number of crystals with varying thicknesses. 

The voluo -of A of KKbO^ determined on the diffrac- 
tometer is 55.934 cm"'’, in fair agreement with the theoretically 
oaloulatcd value of 58.207 om-i Ihe value of/iolserved from 
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the diffractometer is used in the calculations, since (i) dif- 
fractometer gives a "better agreement with the calculated value 
than the photographic one, and (ii) the accuracy of intensity 
measurement is 1% in the diffractometer method compared to 
aioout 3% iu the camera techniq_ue. 
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iPPMDIX IV 


De^bye ten^jerature of Potassium liobate single crystal 

Debye temperature is one of the most important physical 
quantities of a solid, for example, it gives a measure lAhere 
quantim corrections begin to become important for each substance 
and provides useful information regarding vibrational spectrum 
in a qualitative way. Very recently Sastry et al (1969) des- 
cribed how one can find out the energy of formation of defects 
(vacancy) in a solid if the Debye temperature of ''substance is 
known; besides this, Debye temperature is also a very useful 
parameter used in some thermodynamic relations. Potassium Hio- 
bato belongs to an orthorhombic crystal class ^mm2i , and ex- 
plicit formulae are obtained for the Debye characteristic tenpe- 
rature, by an extension of Houston’s method (1948) described in 
detail by Joshi (1961). 

The Debye theory of specific heat is based on a continuum 
model i.e. where the wavelength of the elastic waves is equal 
to or larger than the intcrrtomic distances. This condition 
is satisfied strictly only at the low frequency region of the 
vibrational frequency spectrum. The low frequency waves are 
most important at temperatures approaching 0°K. One would 
therefore, expect that the Debye temperature (9^) derived from 
specific heat measurements should agree with the Bjj obtained 
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from elastic constants, since the wavelengths used for the deter- 
mination of the elastic constants are large compared to the inter- 
atomic distance. 


Many methods of evaluating 9 (elastic) for cubic (Blackman 
1955, Alers and Neighbours, 1959), hexagonal (Wolcott 1959), tetra- 
gonal and trigonal (Betts, Bhatia and Horton 1956) are described. 

An analogous method for orthorhombic crystals is given by Joshi 
(1961). 


In Dehye continuum model 9 (elastic) is given by 


e 


elast ic 




9H N 


1/2 


(lV-1) 


where N is the number of vibrating units in the volume ?, h is 
the Planck's constant, k = Boltzmann constant and I is defined by 

(lV-2) 


3 

2 / 
i=1 0 


lit ^ 

3 


where, v^ axe the Telocities of propagation of elastic waves in 
the continuum at 0°K and dJl-is the solid angle for the whole 
sphere. These velocities are given by three roots of the third 
order Christoffel equation (love 194-4) and hence, in general, 
involve the solution of a cubic equation. These roots will be 
functions of elastic constants and direction of propagation of 
the wave. The problem of evaluation of 9(elastic) thus reduces 
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to the evaluation of the above integral I. 


Si. IV-2 becomes, 


I = ^ (lY-5) 

^ 0 
•3 - 5 /2 

since f(0, 0) = Ej_| Vj_ j ^ , where is the density of the 

crystal. f(9, 0 ) may be expanded in terms of harmonics, P-i _( 0 , 0 ) 

J. I xo. 

having the same symmetry as the Corresponding Ohristoffel eq_ua- 
tion (Betts, Bhatia and Horton 1956), 


OO X 

f(e,0) = s E 

1=0 m=-l 


^m ^Im 


(lY-4) 


Stom the properties of spherical harmonics we have the exact 
relation, 

/ f(e, 0) diL = a^Q = Sq (IY-5) 

O.ur approximation consists in stopping the summation in Bl. IY-4 
at a particular 1 and m. Ihe coefficient a^ is found by solving 
the linear equations IY-4 obtained by giving various values to 
( 0 , 0 ). 


For orthorhombic crystals we take a two fold' rotation axis 
as the polar axis, 9=0. The sGhristoffel equation is invariant 
under the transformations, 

e = Tc-e, 0 = -0 or 0 = 0+71 . 


The appropriate ‘orthorhombic harmonics' are 

°21,2in(®' ^ 


(lY-6) 



188 


With _< 21, 1, m = 0, 1, 2, solving the corresponding equa- 
tion (lV-4), we get, 


% = / "^5 (lY-7 i) 

Uq = [I20f^+21f3-48f^+48f^-f64%+75%l / 180 ClY-7ii) 

^0 " C519f_^+45f3-1008fQ+480f35+1344%+1575%+825%]/3780 


(IV-7iii) 

in vAich the subscripts on f signify its value in the correspond- 
ing directions. The directions chosen are A(100), B(00l), C(010), 
D(101), E(1T30), E(021) and G(102). In deriving (IV-7i) , (IY-7ii) 
and (lY-7iii) all orthorhombic harmonics (IY-6) upto 0^2»®44 
OgQ respectively, are utilized, since at O^q equations (lY~7i), 
(lV-.7ii)and (IV-7iii) are found to he sufficiently convergent. 

The Chistoffel equation for the sound velocities for general direc- 
tion (1, m, n) in the case of orthorhombic crystals, takes the 
form, 

Ciil^-+Cj.gm^+G55n^- <^v^ {C^^+Q^^)lm 

■* § 

(0i2+066)lm Gggl^+Cggmho^^nh (023+0^^)® 

(0^3+055)111 (C23+0„)im 055l2+c^/+0332i2- 

On solving this equation along the directions A, B, ...G, defined 
above, one gets the following expressions for f^^, fg etc. 


-3/2 ^ ^-3/2 

c n + n C 55 3 


-,- 3/2 

+ 11^66 3 

P -,-3/2 

33 -i ' j- - 44-1 L^55-‘ 




contd, . 
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+033+2055)+! 

+E 5 :(Oii+ 033 + 20 j 5 )-|!^(c^^-Ojj)^+ 4 .(a^j+Ojj)^j''/^" 5/2 

. *.3 J 2 

3044+055)3 +l 5 ( 0 ^ 1 + 3 Cg 2 + 40 gg)+'^|(C^^- 3 ag 2 + 2 Cgg)® 

+1 2 (C^2+®66 ) +3022+^66 ) 

- S {(011-3 022+2066 )^+12(0i2'^°66 

Ip=Ir 5 '( 055 +W 66 )] ^'^^+Bfe('^ 022 + 033 + 5 C^^)+:i {(033+3044-4022)^ 

+16(023+044)031/5-3/2 +£-^(4022+033+5044) . 

-1^ {( O33+3G44-+C22) 2+16 ( 033+044) 0 jl/^ - 5/2 

VE 5 (^O 44 +Ogg)J- 5 / 0 +X:^rOii+ 4 O 33 + 5 O 35 )+Tl {(0^1+3035-4033)0 
+16 ( 0i 3+C53) 0 jl/2 5 -3/24X-^ ( 0i 1 +4C33+5C55 ) 

-'iS {(Oii+3O55-4O33)O+l6(Oi3+O33)0jl/05-3/2 (IT-9) 

The characteris'tiG temperature for orthorhomhic Potassium 
Kiohate can he calculated from the formula (lY- 7 iii). Eq.s. (IY- 7 ii) 
and Cl¥-7iii) differ only little from each other (< 2%),(Joshi, 
1961 ). One has to use higher term approximations than those used 
here to improve the accuracy. 

Thr an accurate evaluation of 0( elastic), the values of 
elastic constants measured near 0°K should have been used. In 
the absence of such measurements, elastic constants measured at 

he used, but, then, the absolute values 


other temperatures may 
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of the characteristic temperature becomes xmreliable to some 
extent. Howe-ver, Blackman (1955) has pointed out that the varia- 
tion in 9 (elastic) values is not very large relative to the varia- 
tion of elastic const arts with temperature. Uo specific heat 
data are available for Potassium Niobate in the true region to 
estimate 6^. 

0 (elastic) for Potassium liobate = 302°K 

This value is obtained from (lV-7iii) which is expanded 
upto 6o terms i.e. until it is fairly convergent. 

Inversion Temperature: 

The Bragg intensity (corresponding to reflection from 
ideal planes) decreases with increasing T because thermal vib- 
rations make the plane less and less perfect. Diffuse scatter- 
ing increases with temperature since thermal disorder increases, 
but it increases only up to the point where the plane retains 
its physical meaning relative to the effect being studied (i.e., 
first, second and third order diffuse scattering). Thermal 
diffuse scattering (TDS) begins to decrease i/dien, as & result 
of thermal agitation, the plane begins to lose its reality, 
that is, the crystal starts to 'melt'. The teniperature where 
this takes place is ealled the inversion temperature, later 
the whole structure melts, and the crystal ceases to exist. One 
can find out the inversion temperature (l^ax^ order 
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diffu-e scattering (TDS-^) for a given plane, if the Debye tem- 
perature of the crystal is known, from the following relation 
(Amoros 1968) 

= -f ? ”4 4i . 

Symbols have their usual meaning. 


The inversion phenomenon is quite general and the ex- 
pression for this is independent of the type of structure, and 
is a function only of the mean square amplitude of atomic vibra- 
tion, u^ which ia in turn related to the interplanar spacings, 
(dy^j^l), therefore only certain sets of crystal planes will be 
potentially able to show the inversion phenomenon. 


Relation between atomic- vibrational amplitude and inversion 


temperature) : 

The atoms in the crystal vibrate with characteristic 
amplitude for a given structure and temperature. Since the phe- 
nomenon of diffuse scatterihg is a direct consequence of this 
atomic vibration, there must exist a relation between the atomic- 
vibrational amplitude and ' 

Letting T = (T^^^) ^ , one can obtain from the well loiown 
Debye- Waller relation (James 1954, imoros 1968)^ 


(2 


max^ 1 „ 
T 





0, 16 


(IV-11) 


and 


(IV- 12) 
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Therefore, TDS^ (1st order thermal diffuse scattering) 
reaches its maximum value for a given crystallographic plane 
when (u ) is a definite fraction 0.16 of its spacing, Grune is en 
(1925) showed that melting takes place when (u^)^'^^ is of the 
order of 0.10 of the shortest interatomic distance. These 
relations are simple to derive for a cubic monoatomic crystal. 

But for the non cubic polyatomic crystals, these relations can 
be obtained by considering the mean square amplitude of vibration 
of each kind of atom separately, as well as the anisotropy of 
its vibration. 
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EREATA 




O 


G3 


1Q> 


B1 


81 


Line 

linos 16 and 18 

5in line from Lot tom 
(Reference) 

14 

lines 1 and 2 


Fig. 4.6(a,'b) 
(Caption) 


Correction 

" Curve corresponding to” instead 
of ” Curve corresponds to” 

” J. PLys. Rad. ” instead of 
” J. Rad. PLys. ” 

” CLakraLorty ” instead of 
” Sen and CLakraToorty” 

” adjusted at 0.02 cm in LeigTat” 
instead of ” adjusted at 0.02 cm 
to 0.02 cm in height” 

” The curve corresponding to” 
instead of ” The curve corres- 
ponds to” 


Fig. 4.6(a) 
(Caption) 


insert "(as shown in Fig. 4_,6h)” 
before ” is ” 

2 

” of cross section 5 x 5 nim ” 


106 


21 


instead of ” of the size 
5 xm X 5 nm” 





